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5^ ' We set constraints on moduli cosmology from the production of dark matter - radiation and 

Ci^' baryon - radiation isocurvature fluctuations through modulus decay, assuming the modulus remains 

light during inflation. We flnd that the moduli problem becomes worse at the perturbative level 
as a signiflcant part of the parameter space (modulus mass) - tJinf (modulus vev at the end of 
00 ' inflation) is constrained by the non-observation of significant isocurvature fluctuations. We discuss 

in detail the evolution of the modulus vev and perturbations, in particular the consequences of 
Hubble scale corrections to the modulus potential and the stochastic motion of the modulus during 
inflation. We show, in particular, that a high modulus mass scale 771^ > 100 TeV, which allows the 
modulus to evade big-bang nucleosynthesis constraints is strongly constrained at the perturbative 
. level. We find that generically, solving the moduli problem requires the infiationary scale to be much 

^' ■ smaller than lO" GeV. 

o 

^\ 

I. INTRODUCTION 

(N 

^ ' The existence of scalar fields with gravitational coupling to the visible sector appears to be a generic prediction of 
\^ \ particle physics beyond the standard model. This, however, may cause serious cosmological difficulties, as exemplified 
CN ■ by the "cosmological moduli problem" [l|, Qj- Assuming that the mass rua of a modulus is of order of the weak 
' scale, as one would expect for soft masses induced by supersymmetry breaking, this field should decay after big- 
^— ; , bang nucleosynthesis, on a timescale ta ~ 10 sec (TO^/100GeV)-3/2. The ensuing high energy electromagnetic and 
' hadronic cascades would then ruin the success of big-bang nucleosynthesis predictions (see 0] for a recent compilation 
, and references therein) unless the modulus energy density were extremely small at that time. By "extremely small" , 
it is meant about twenty orders of magnitude smaller than what is generically expected for a scalar field oscillating in 
a quadratic potential with an initial expectation value of order of the Planck scale. Turned around, this cosmological 
moduli problem reveals the power of big-bang nucleosynthesis when used as a probe of high energy physics and early 
Universe cosmology. 

^ Several classes of solutions have been proposed. The first one argues that the vacuum expectation value (vev) 

d ' CTinf of the modulus at the end of inflation is much smaller than the Planck scale 0, [E B S 111 ■ This is not a 
trivial requirement as it demands that the effective minima of the modulus potential at low energy (i.e. well after 
post-inflationary reheating) and at high energy (i.e. during inflation) coincide with each other 0, Q. Furthermore, 
quantum fluctuations of the scalar field will generally push the field away from this minimum Q . 

An alternative solution to the moduli problem proposes that the modulus mass is so large that the modulus decays 
before big-bang nucleosynthesis, leaving enough time for the high energy cascade to thcrmalize before the process of 
nucleosynthesis actually starts [lOl [Til. [Tsl Il3| . This requires mg. > 100 TeV. Although this lies some two orders of 
magnitude beyond the expected soft scale, such masses can be accommodated in successful models of supersymmetry 
breaking such as anomaly mediation or no-scale supergravity as argued in the above references. 
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Finally, it has also been proposed to dilute the energy density contained in the moduli through an epoch of low 
scale inflation [l^ or thermal inflation [l^, [l3| • 

As formulated above, this standard moduli problem is directly expressed as a constraint on the energy density 
of the modulus field at the time of its decay. Meanwhile, progress in observational cosmology has been such that 
it is now possible to constrain the nature of density perturbations to a high degree of accuracy. Most notably, the 
analysis of microwave ba ckg round t emp erature fluctuations allows to constrain the fraction of isocurvature modes to 
a quite low level [H, [H, Ho, [U, 111, Hi [H, HE M, [13, HI HI, S HH . As we argue in this paper (see also ) , such 
constraints on the spectrum and the nature of density perturbations can be translated into constraints on moduli 
cosmology. We will find that the moduli problem becomes worse at the perturbative level. The main reason is that a 
modulus, being uncoupled to fields of the visible sector, inherits its own fluctuations through inflation. At the time 
of reheating, there exists an isocurvature fluctuation between the modulus and radiation, which is transformed into a 
dark matter - radiation or baryon -radiation isocurvature mode when the modulus decays into radiation, dark matter 
and baryons. In this way, one may thus picture the modulus as a curvaton field, whose phenomenology has been 
intensively scrutinized in the past few years (see notably Refs. [sl, H S 111, [S^, [SI, 111, SO, [S S 

The generation of isocurvature fiuctuations by modulus decay has been noted before, see ST], but to our knowledge, 
a detailed analysis of the ensuing constraints on moduli cosmology has not been given up to now. The present paper 
furthermore attempts at being general and exhaustive with respect to modulus cosmology. In the course of our 
discussion, we have thus obtained new results in several places, such as those related to the evolution of the modulus 
and its perturbations in the presence of supergravity corrections to the modulus potential, or generalized existing 
discussions, for instance concerning the stochastic behavior of the modulus during infiation. The layout of the paper 
is as follows. In Section |TI1 we describe in general terms how constraints on isocurvature fiuctuations can be turned 
into constraints on moduli cosmology. In Section [IIII we then calculate for various modulus effective potentials (time 
independent, or accounting for supergravity corrections) the cosmological consequences and present the constraints in 
the modulus parameter space (mo-, a-mi)- We summarize our findings and discuss how to evade the modulus problem 
at the perturbative level in Section |lVl Finally, the paper ends with three technical appendices, which contain results 
of importance to the present study but that also possess interest of their own. The first one. Appendix \X[ is devoted 
to the calculation of the quantum behavior of the modulus field during infiation, for large field and small field models, 
accounting for possible supergravity corrections to the modulus potential. The second one. Appendix iBl discusses in 
greater details the evolution of the modulus vev when its potential receives supergravity corrections. The third one, 
Appendix[Cl presents some supergravity based, concrete particle physics models for the modulus field. All throughout 
this paper, Mp\ — 2.42 x 10^^ GeV denotes the reduced Planck mass. 



II. GENERATION OF ISOCURVATURE PERTURBATIONS 



This section sets the stage for the next Section, which performs a systematic study of the constraints obtained in 
moduli parameter space. Here, we introduce the relevant physical parameters and we describe how the amount of 
isocurvature fluctuations produced through modulus decay can be calculated analytically. The formulae obtained will 
be useful to interpret the results of numerical calculations presented in the following Section. 



A. Background evolution 



Let us first sketch the cosmological scenario and the outline of the calculation. We assume that infiation proceeds at 
the energy scale Hini- At the end of infiation, the slow-roll conditions are violated (or, in the case of hybrid inflation, 
an instability occurs), the inflaton field 4> starts oscillating rapidly at the bottom of its potential and the reheating 
period begins. If the infiaton potential is quadratic, then the universe becomes matter-dominated. At a later stage, 
the inflaton field decays into radiation, dark matter and baryons. At this high energy scale, well above the dark matter 
mass, baryons, radiation and dark matter are all part of the same "radiation" fiuid and the universe is effectively 
radiation-dominated. The temperature T^h of the radiation fluid at the beginning of this era, the post-inflationary 
reheating temperature, is a direct function of the decay rate of the inflaton field. 

Let us now consider the modulus field a. In the following, we denote its vev at the end of infiation by Cinf and 
treat this quantity as a free parameter (however, see below the considerations on the quantum behavior of a during 
infiation). In the post-infiationary era, in order to follow the evolution of cr, the shape of the potential is needed. 
In what follows, one considers two cases: one where the potential is purely quadratic and one where the potential is 
affected by Hubble scale contributions, meaning that a term of the form [a ~ ao)^ is added to the quadratic part. 
If the potential is purely quadratic and if the modulus field is a test field, then its vev remains constant (therefore 
equal to Cinf) until H = nia- Let us notice that, when H = ma, the energy density of the background is SMp^m^ 
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while p„ ~ "2^(T?jf. Then, the condition that cr is a test field, that is to say that its backreaction on the expansion 
rate is negligible, i.e. < p^, implies that Cinf < Mpi. We thus do not discuss the possibility raised by Linde and 
Mukhanov Q that the modulus becomes itself an inflaton after having acquired a large vev through quantum jumps. 
On the other hand, if the Hubble scale corrections are important, a can never be considered as light (unless c <C 1) 
since its effective mass is always of the order of the Hubble scale. As a result, the evolution of the modulus between the 
end of inflation and the time H — nia can become rather involved since a has no reason to stay constant anymore as 
it was the case in the purely quadratic situation. For this reason, we postpone the detailed discussion of the modulus 
evolution to the following section and we always express our results in terms of the modulus energy density when 
H = TOct- 

At late times H < m^, Hubble corrections to the mass term indeed become, by definition, negligible. As a 
consequence, the modulus potential is then given by At iJ = m^, the modulus starts oscillating at the 

bottom of this potential and pa- cx a^^. This occurs at an equivalent temperature scale given by 

r.., = (il!|p) , io"GeV . (1) 

Of course, Tosci can correspond to the temperature of the radiation bath only if Tosci ^ T^h- Otherwise, it should 
be thought of as the temperature that the radiation bath would have were the energy density contained in radiation. 
Typically, one has Tosci > Tih unless the reheating temperature is very high. At the onset of oscillations, the a field 
carries a fraction r^o^.osci of the total energy density which can be expressed in terms of (Tinf , the vacuum expectation 
value of a at the end of infiation. 

Modulus eventually decays into radiation, dark matter and baryons when H = T^, T^r being the gravitational decay 
width of cr: 

As a consequence, the decay of the a field occurs at a temperature: 

- (r.„.,,- . lO-GcV {t^y'"{^f . ,3, 

Clearly, the decay occurs much after the onset of oscillations and the reheating. 

Finally, let us end this short description of the background evolution by mentioning that we assume all throughout 
this paper that dark matter originates from freeze-out of annihilations. The dark matter freeze-out occurs at a 
temperature of m^/x[, where is the dark matter particle mass and X[ 20 — 30. Since, typically, 
O (100) GeV, one obtains a temperature of ~ 1 — 10 GeV. Therefore, provided rricr < 10^ GeV, the modulus always 
decays after dark matter freeze out. On the other hand, the freeze-out of the baryons takes place at a temperature of 
~ 20MeV. Therefore, whether the modulus decay occurs before or after the baryons freeze-out depends on the value 

of JTicr. 

The previous considerations imply that two crucial variables in this study are and ainf and, therefore, in the 
following, we will express our constraints in the plan (mo-, cri„f). The mass scale should in principle be fixed by 
high energy physics with a preferred range around 10^ — 10^ GeV. On the contrary, the vev Cinf is determined by the 
early cosmological evolution. The mass scale determines, among others, the decay time of the modulus, and together 
with the vev (Tinf, it also determines the magnitude of the modulus energy density at the time of decay, hence the 
amount of isocurvature perturbations transfered to the dark matter and baryon fluids. The modulus vev ai^f can 
be expressed as the sum of two parts, one corresponding to the classical trajectory ad of the modulus field in its 
potential during inflation, and the typical spread (Ja^)^/^ around this trajectory due to quantum effects [isl. liol. [sot . 
The standard deviation on scales larger than the Hubble radius (Sa^)^^^ has been discussed recently by Linde and 
Mukhanov [9|] and Lyth [51]. It is discussed in greater detail in Appendix [XI along with the classical trajectory a^i 
during inflation. 

For both large field and small field models, one may summarize the situation as follows (the case of hybrid infiation 
is also treated in Appendix [XJ. Consider first the case in which the modulus is effectively massless during inflation, 
meaning: 

-^inf,in / ,\ 

ma,eS < , (4) 

where i?inf,in is the Hubble parameter at the beginning of inflation and iV^ the total number of e-folds of inflation. 
Let us notice that, strictly speaking, this condition is not equivalent to TOo-.cfi ^ -ffinf if the Hubble parameter evolves 
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during inflation and/or if the total number of e-folds is large. The modulus mass is written mo-. off to encompass two 
different cases: a fixed mass m„ or a (supergravity induced) Hubble scale mass CjT?, with < Cj < 1. Here, we have 
written Ci in order to emphasize (see also Appendix [A] and see Appendix [C] for concrete examples) that the effective 
mass during inflation is not necessarily the same as the effective mass in the post-inflationary era (in other words, a 
priori, Ci ^ c). For large field m'^cfP' infiation, ifinf,in-/V~^^^ = \/ll'im^ ~ 1.4 x lO^'^ GeV is fixed by normalization to 
the cosmic microwave background anisotropies. Given the values < 10^ GeV considered in this paper, this means 
that, in the pure quadratic case, the modulus field is always massless. Of course, when Hubble scale corrections are 
present and Cj ^ 1, this is no longer the case. For small field infiation, however, the situation is different. Indeed, 
the quantity i?inf,inA^^^^^ can be very small or very large depending on the infiationary scale (which can be as low 
as ~TeV) and the number of e-folds {N^ > 60 but is otherwise essentially unbounded, see Appendix E]) . Therefore, 
for small field infiation, the pure quadratic case can or cannot correspond to a massless situation. At the classical 
level, one finds that, in this "massless" field case, CTcI ~ crin, meaning that the classical value has not changed during 
inflation. 

Concerning the contribution of quantum effects, one finds: 

2N1/2 ^ :HinUn^l/2 (5) 

This result holds for small field infiation; for large field to^0^, it is a factor \/2 smaller, see Appendix [XI This value 
does not depend on ma,eS and it diverges in the limit of de Sitter spacetime {N^ +00), as expected for a massless 
field. This value can actually be understood simply as follows: every e-fold Hinft, the field performs a random step of 
length ±7?inf / (27r) , which add up randomly, yielding the above random walk behavior. Setting > 60 yields the 
following lower bound: 

(^.y^>5xlO-Mp,(^||g^) . (6) 

For small field inflation, i?inf,in ^ Hmf, i.e. the Hubble constant does not change much during infiation. For large 
field rn^cj)^ infiation however, as already noticed before, the numerical prefactor is 1/a/2 times smaller, but i?inf,in — 
y/2N^Hinf, and Hinf ^ 10^^ GeV, so that overall the above bound is a factor ~ 8 larger, see Eq. (|A33p . 

Let us now consider the other limit of a massive field, mo-. eft 3> Hi^itsnN^^^'^ , yet not too massive in the sense that 
Wa-,cff ^ Hinf. This is relevant for the pure quadratic case during small fleld inflation (depending on the parameters, 
see above) and when the Hubble scale corrections are present in both cases. Then, at the classical level, one finds 
that the fleld evolves and rolls down its potential during inflation. Regarding the magnitude of quantum effects, one 
obtains: 



{SaY' ^J±^JS^. (7) 

For large field m'^clP inflation and typical moduli masses ^ 10^^ GeV, this case only applies if the modulus receives 

Hubble scale mass corrections mcr,cff = CiH, since Hin{,inN~-'-/^ ^ 10^^ GeV. Furthermore, the right hand side in 

the above equation should be multipHed by an extra factor of ^/3 + c? /q in this case, see Eq. (jA44[) . Of course, 

as Ci ^ the fleld becomes light and one recovers the previous result, see Appendix El In the case of small 

field inflation, the above value can be quite small and all the more so as the scale of inflation is lowered. Note 

1/2 

that the above value of {dcr'^') reproduces the well-known Bunch-Davies expression for a massive field in de Sitter 
spacetime [s^, [H, [U, [S^ . As discussed in Ref. jsl , one can understand this result by considering the same random walk 

as before, but noting that modes on large wavelengths redshift away in proportion to exp — ^jji/(3-ffinf ) , which 
implies that the maximum contribution to the fluctuations has been generated during the last AN ~ 3H?^f/m^ 
e-folds. The product AA^ x H?^f/{4'K^) then reproduces the Bunch-Davies result (squared) to within a factor 2. In 

this respect, one should note that the previous limit mo-. oft ^ Hi^i / \/N^ that we considered corresponds to a field so 
light that modes do not have time to redshift away in N,j, e-folds. In this limit, AA^ is bounded by A'^, hence Eq. ([5]) 
is recovered. 

Finally, the last case of interest is mo, eft 2> -ffinf. In this situation, the fleld is too massive to be excited, and 
1/2 

consequently (^cr^) is exponentially suppressed. 
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B. Evolution of perturbations 



Let us now introduce the scenario at the level of perturbations. As we have just done for the background quantities, 
one can also follow the perturbations of each species throughout the cosmic evolution. To be more precise, we are 
interested in the curvature perturbation for the species "a" defined by [s^, [13, HI, 



C^^-^-H^c.-^+ (8) 

Pa 3(1+Wq) 

where $ is the Bardeen potential, Aq, the gauge-invariant density contrast and uja = Pa/ Pa the equation of state 
parameter. 

After the decay of the inflaton field, the fluctuations in (j) have been transmitted to radiation, characterized by , 
dark matter, C,^^ and baryons (and anti-baryons), Cg''- Since these fluids share thermal equilibrium, one has [63| 

C}y^ = C,^ = c}^^ = C,^\ Let us notice that, at this stage, dark matter and baryons are still relativistic fluids. Indeed, 
dark matter becomes non-relativistic at a temperature of ^ = O(IOO) GeV. Regarding the baryons, the situation 
is more complicated since, in principle, they become non-relativistic at a temperature of ~ 1 GeV (about the same 
as the dark matter freeze-out temperature), that is to say well below the reheating temperature, but in fact at that 
temperature one still have a quark-gluon plasma. 

We define the above initial conditions, indexed with well into the modulus oscillations era, at <C rria and 
T < Tj-h, and before the modulus comes to dominate the energy density. In this era, the modulus can be considered 
as a pressureless fluid, supergravity contributions to its potential have become negligible, hence previous results 
on curvaton phenomenology can be applied, as discussed further below. One needs however to relate the modulus 
curvature perturbations at this time, C,a^ to the modulus perturbations acquired through inflation. This obviously 
depends on the modulus potential at rria < H < Hinf- 

For a simple time independent quadratic modulus potential, one can use the results of Ref. which give: 

C«=-|*in.+ |^. (9) 
Z 6 (Tinf 

The quantity (E>inf denotes the Bardeen potential at the end of inflation, Sainf denotes the modulus perturbations on 
large scales and the calculation assumes that the modulus behaves as a test field, so that $ is approximately constant: 
^inf — It also assumes that radiation dominates at the time at which is defined. It is important to realize 
that the radiation curvature perturbation is related to the Bardeen potential through C^'' = — S^^'V^, so that the 
initial modulus - radiation isocurvature perturbation can be rewritten as 

5« - 2^ . (10) 

Cinf 

In Section IIIIBl we show that this result holds even when the modulus potential receives a supergravity inspired 
+c^H^ quadratic mass term. This result is of importance for the present discussion, since it shows that the modulus 
- radiation isocurvature fluctuation disappears in the limit ^Cinf/cinf — *■ 0. One way to achieve this is to assume that 
the modulus is heavy during inflation, either because iJinf ^ rn^ or because the modulus receives an effective mass 
term +cfH^ during inflation. Furthermore, Ref. [6l| has shown that the isocurvature mode between the inflation and 
any heavy field actually disappears during inflation because the heavy field is drawn to the minimum of its potential 
at every point in space, so that in this case, there would not even be an isocurvature fluctuation to start with, at the 
end of inflation. For this reason, we discard for now this case and assume everywhere that the modulus has remained 
light during inflation, in which case ^(Tinf/cinf — -ffinf/(27rcrinf). In Appendix [Cl we present several concrete models 
of inflationary model building in a supergravity framework; for both models of D-term inflation, it is found that the 
modulus remains light during inflation, but acquires a Hubble effective mass after inflation. 

In the following, we set C-y''' — 10~^. From the time at which the initial conditions are defined, all the (a remain 
constant until dark matter freeze-out (in between, the dark matter has become non-relativistic, see before). During 
this phase, the radiation, baryons and modulus fluid perturbations are not affected, 

c>' = c^^ e- = c^'\ c^' = c^. e = c«, (11) 

but the dark matter perturbations are modified according to [4 

(«f - 



C = c 



.(i) _ ^(i) 



(12) 
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with ai = xi + 3/2. From the above equation, one sees that the quantity Cx is not modified if (i.e. the 



^(i) 



in which case the freeze-out surface exactly coincides 



modulus is negligible at dark matter freeze-out) and/or (a^^ = 
with the uniform radiation surface. 

Then, the ^q's remain constant until the baryons freeze out (assuming it occurs before modulus decay). Through 
this stage, it is clear that C^, and Co- remain unaffected. On the contrary, one expects (^b and to evolve. One can 
consider a "net baryon number" fluid, the energy density of which is given by f^b — ^S- Before baryons freeze out, this 
fluid is made of baryons and anti-baryons in thermal equilibrium (with a small excess of baryons) but after the freeze 
out of annihilations, it is essentially made of baryons. In the absence of any baryon number violating process, this 
fluid of "net baryon number" is isolated, hence its curvature perturbation remains constant. Therefore, after baryons 
freeze-out, one has Cb = Cb"*^ — Cb and Cj, = 0. Notice that the same reasoning would also be valid in the case where 
the freeze-out occurred after modulus decay. In fact, the above discussion would be modified only if baryon number 
violation occurred after modulus decay. 

Finally, the modulus decays in dark matter, radiation, baryons and anti-baryons and it is clear that all the corre- 
sponding curvature perturbations are then modified. We obtain (see Refs. [45„i47j for details): 



= c<^ + 



/■<d _ ^(i) 



7(C 



<d 



^x J 



So- Sy 



(13) 
(14) 
(15) 



As explained in Refs. [H, [43|, in order to obtain these relations, we have assumed that a fraction B^ = T^^/Ta- ^ 1 
of the cr energy density goes into dark matter particles. One can relate this parameter B^ to the branching ratio of 
modulus decay into dark matter particles as follows. Assuming that dark matter particles thermalize instantaneously, 
then if one modulus produces through its decay particles, one finds that B^ = N^m^/m^. In the range of 
parameters that we are interested in, <C nia and < 1, see Refs. [H, |6l| for a detailed discussion. Strictly 
speaking, the thermalization is quasi-instantaneous only in the high modulus mass range, ma ^ 100 TeV, while at 
lower masses some rcdshifting due to the cosmic expansion occurs. This can be seen as follows. Using the results of 
Refs. [H, nil, one can write down the scattering cross-section of Compton-like processes x + ^ ^ X + ^ (with i an 
ultra-relativistic lepton of the thermal bath, e.g. a neutrino or an anti- neutrino at the BBN epoch): 



3C {i 



1287r 



9 1 



(16) 



with C a prefactor of order unity deflned in Ref. [64| . and s the standard center of mass energy squared: s ~ 
2E^Ei ~\- with ~ 3.15r where T is the temperature of the thermal bath. In the region of interest (at modulus 
decay), one can check that 2E^Ei ^ m^. The ratio of the interaction rate, T^^i = niaxi to the Hubble rate at the 
time of modulus decay can then be expressed as 



H 



1.4 X 10*^(7 



E^ 



100 TeV 



13/2 



100 GcV 



(17) 



As a consequence, if ma ^ 30TeVC~^/^ (m^/lOO GeV)^^^, one has T^e > H for all energies i?^ > m^, which implies 
that the x-particle becomes non-relativistic through multiple interactions in less than a Hubble time. In this range, and 
for our purpose, one can treat the dark matter fluid as a pressureless fluid immediately after modulus decay. Inversely, 
if ma ^ TTeVC"^''"'^^ (m^/100 GeV)^^''"'^^, F^^ < H for all energies E^ < ma/2: the particle never interacts and simply 
redshifts to non-relativistic velocities within \ii[ma / {2mx)] e-folds. Finally, in the intermediate range, the x particle 

looses its energy through interactions to some intermediate value 90GeYC-^^^ (m^/lOO TeV)"^/^ {m^/lOOGeVf, 
then redshifts away down to m^ through cosmic expansion. 

All in all, our above assumption of "instantaneous thermalization" amounts to neglecting this redshifting factor, 
which in turn slightly overestimates the abundance of a produced dark matter by a factor which never exceeds ~ 35. 
This value applies at ma — 7 TeV and it rapidly decreases to unity away from this value. This only affects very 
marginally the results derived below. 

We have also assumed that the fraction of modulus energy density Sb+b that goes into baryons and anti-baryons 
is very much smaller than unity, as one would expect. With respect to Ref. [47[, we have also assumed here that 
modulus decay preserves baryon number. Finally, the parameter r*^'' that appears in the first of the above equations 
has been introduced in Ref. [36l|; if dark matter is entirely produced by modulus decay (i.e. fl^' ~ 0), then 1 — r*-"* 
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characterizes the amount of initial modulus - radiation isocurvature mode that is transfered through modulus decay. 
We use the simple formula: 

which has been found numerically to be a good approximation f4ll| . 

After modulus decay, the Cq's remain constant throughout the subsequent cosmic evolution. The corresponding 
values can be compared to CMB data. 



C. Transfer of isocurvature modes 



As we will see shortly, dark matter - radiation and baryon - radiation isocurvature modes are generated in different 
regions of the parameter space. The constraints obtained are thus complementary to each other. For this reason, 
and for the sake of clarity, we discuss the generation of each mode in turn. Let us also recall that the isocurvature 
perturbations between two fluids "a" and is defined by Sa/3 = 2> {C,a — Qp) ■ 

The dark matter - radiation and baryon isocurvature modes on large scales are given by [45j |: 



C>d 



1 



1 + T 



X 



s^l , (19) 



c -f}<<^5«. (20) 

These formulae can be straightforwardly deduced from the results quoted in the previous section, except the presence 
of the parameter in S^^ . This parameter represents the ratio of the dark matter annihilation rate to the expansion 
rate immediately after the decay of the modulus field. If this latter produces sufficiently many dark matter particles, 
these may annihilate with each other. As 3> 1, meaning that annihilations are effective, the isocurvature 
perturbation transfer is partially erased (see Ref. [i^ for details). One finds: 

- ^^m/- 0.076 m^Mpi ^-^^ ^^^^'^ VlOOTevJ UoOGevJ VlOoJ ^'^''^ ' ^^^^ 

The above formula neglects the amount of dark matter initially present; this is a good approximation insofar as the 
annihilation rate of these dark matter particles is very much smaller than the Hubble rate after freeze-out in the 
absence of modulus decay, which would lead to ^ 1 hence to a negligible correction to the equation for the 
transfer of isocurvature mode. Note that the last factor involving X{ in the above equation may be quite large, being 
~ 2.2 X 10^ for Xf = 20 and 5.9 x 10^"^ for X[ — 30. Finally, note that the above formula is only approximate (see 
Ref. [45] for details). The constraints that we derive further below are obtained through the numerical integration of 
the full set of equations of motion and are therefore more accurate. 

A significant dark matter - radiation isocurvature mode is generated if both following conditions are satisfied: 

B^n<'' > 17<<' , 1]<'^ < 1 . (22) 

The former condition expresses the fact that the amount of moduli produced dark matter particles exceeds that 
coming from freeze-out of annihilations, while the latter requires that the modulus energy density is not sufficient to 
affect the radiation content. All in all, this means that the modulus perturbation is transferred to the dark matter 
fluid but not to radiation. One could also imagine that the modulus perturbations are transfered to the radiation 
sector but not to the dark matter sector, thereby generating a net isocurvature perturbations. However, this would 
require either blocking the decay of modulus to dark matter, which is unlikely as the dark matter particle is always 
much lighter than the modulus in the parameter space we are interested in, or, having the modulus decay after matter 
- radiation equality, which is forbidden by constraints on cosmic microwave background distortions. 

At this stage, we need to compute explicitly the parameters appearing in Eqs. (|19p and (|20p . In particular, from 
the expression of the dark matter annihilation cross-section, one obtains 

n>' ~ 1.67 X lO-^Xf /^e-^f . (23) 

One also needs to evaluate the quantity ilj^'' and il^''. From the fact that the energy density of the modulus scales 
as a~^, one obtains that fla^ just after the dark matter freeze out can be expressed as 

= r!,,osci— min (Tosci , T,^) . (24) 
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In this expression, the minimum of Tosci and Tj-h appears because, if the oscillations start before the end of reheating, 
then the inflaton and modulus energy densities have the same scaling until the reheating stage is completed. Only 
below Tosci or Tih, whichever is smaller, the energy density of the a oscillations increases relatively to radiation energy 
density. In order to obtain the above formula, we have also assumed that the modulus can never start a new phase of 
inflation. In the same manner, immediately prior to decay, the ratio of the energy density contained in a oscillations 
to that contained in radiation reads: 

We are now in a position where one can calculate the transfer of isocurvature perturbation from modulus - radiation 
to dark matter - radiation and baryon - radiation. This is done in the Section that follows. 

III. CONSTRAINTS IN THE MODULUS PARAMETER SPACE 

Our present goal is to compute the amount of isocurvature perturbation produced through the differential decay of 
the modulus into dark matter and radiation, as well as between the baryon and the radiation fluid, which can then 
be compared to existing bounds obtained through the analysis of cosmic microwave background fluctuations [H, [l^, 
[20l . I21I . [2^ [23l . [23 . [25I I26L [27l . [28L I29L [30I I . We have chosen to express the amount of isocurvature fluctuation in this 
matter sector as follows: 



Cm C7 
(Cm + C7)/2' 



S^, ^ 7T^^^^^^ (26) 



where the subscript "m" comprises all of non-relativistic matter, i.e. dark matter and baryons (so that, for instance, 
Jim = + f^b)- The quantity Cm can be written in terms of the baryon and dark matter curvature perturbations Cb 
and C^: 



Cm ^ §^Cb + §^Cx . (27) 



The definition (pS)) can be justified by the fact that the data are in fact sensitive to the quantity defined by S^^ = 
Sh-y + i^x'^xi/^b, see Ref. [i^. In this reference, the quantity B = Sg^/C-y is constrained using various cosmic 
microwave background data (including WMAPl). At 95% C. L. it was found that -0.53 < B < 0.43, see Ref. [esj . 
Our quantity Smj is related to B through 



2B 

6nrn7^+ B ' 



•^m^ = ;^ , „ , (28) 



which imphes —0.12 < Sm-y ^ 0.089 where we have taken fl^h'^ ~ 0.12 and Q^h"^ ~ 0.0225. The choice is also 
motivated by the fact that the most recent analysis of Wilkinson Microwave Anisotropy Probe fifth year (WMAP5) 
data has constrained the same quantity for dark matter only (i.e. m — > ;\;) [3l| . These results give an upper bound as 
low as 2.0% (95% C.L.) for fully anti-correlated isocurvature modes and 8.6% (95% C.L.) for uncorrelated modes. In 
the following, we will present contours on the quantity (5ni7 and emphasize the loci of 1% and 10%. 

The quantity ^^-y must be calculated at the time of recombination, well after baryon and dark matter freeze- 
out of annihilations and curvaton decay (and big bang nucleosynthesis). The individual gauge invariant curvature 
perturbations Cb, Cx ^^'^ C7 are then constant since the fluids can be considered as isolated at that time. Given the 
results of the previous section, in order to compute the constraints in modulus parameter space and to evaluate Sm-y, 
the value of rio-,osci, which directly controls the energy density parameter of the modulus at the time of decay, is the 
only remaining quantity which remains to be specified. 

Before embarking on a detailed discussion of these above considerations, one should note that two other cosmological 
constraints are to be satisfied. One concerns the present day abundance of dark matter and the other the overall 
amplitude of the total curvature perturbation. Regarding the former, dark matter is produced both thermally (through 
the freeze-out of annihilations) with present day abundance $7^ q, and non-thermally (through modulus decay) with 
present day abundance fi^ q (immediately after decay fi^ = By,n,^''). Hence, the final abundance is controlled by the 
parameters Xf, m^, and fi^''. In order to minimize the dimensionality of our parameter space, we have chosen 
to proceed as follows. We set X{ — 21 and = 100 GeV, which implies that q = 0.2; we then tune Ny, for each 
value of our main parameters ma and Cinf such that the total ^xfi ^i^s within a factor of three of its observed value. 
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More precisely, we maintain = 1 (see above and Refs. [l^, [g^l) whenever Q'^q < 0.5 and decrease it in order to 
saturate this last bound otherwise. This is somewhat arbitrary, but given the remaining freedom in X[ and m^, one 
could always tune slightly the parameters to achieve a better agreement with the known value. In any case, this 
procedure hardly modifies the constraints obtained in this article, as we have checked. Furthermore, this approach is 
conservative in the sense that we forbid the non-thermal channel to exceed twice the thermal channel which slightly 
reduces the isocurvature perturbations. In this way, in all of the parameter space scanned in the subsequent figures, 
the dark matter abundance is correct to within a factor of two to three. 

Concerning the overall amplitude of the curvature perturbation, one needs to require that ^^'^ ~ 10^^. This 
quantity is determined by Eq. (|13p and it can be rewritten as: 

C>^ ^ C^') + if^<^5« . (29) 

Therefore the magnitude of the total curvature perturbation is controlled by several parameters, including dint and 
Hinf which determine the scalings of ri^"* and S^j. In principle, it would be possible to rescale i?inf in order to reach 
the correct magnitude for (^^"^ at each value of Cinf . However, this would make the interpretation of the figures rather 
complex. In the following, we have rather chosen to plot the constraints obtained for two values of -ffinf in each case, 
in order to gauge the effect of Hinf on these constraints. These two values are Hini = 10^'^ GeV, which provides a 
natural scale for inflation since it corresponds to the simplest inflaton potential m^cj)^ (see also the discussion about 
naturalness in Ref. [66]), and ifmf = 10^ GeV. The latter is chosen arbitrarily, but it is such that the total curvature 
perturbation is of the right order of magnitude at every point of the modulus parameter space. For the former 
value of i?inf , a significant region of modulus parameter space is excluded by the normalization of the total curvature 
perturbation; however, this region is entirely contained in the region which is excluded by the isocurvature constraints. 



A. Quadratic potential 

Here, we assume that the potential of the modulus is a simple quadratic potential V{ij) = rn^a^ /2 from the end of 
inflation onwards. At the onset of oscillations, the a field carries a fraction r2a-,osci of the total energy density: 

Therefore, one can now explicitly evaluate the quantities ri^', fl^"^ ~ r^"*, see Eqs. (P^. (IMl) . (pS)) and In 
particular, using Eqs. ((25)) and ([3]), immediately prior to decay, the ratio of the energy density contained in a 
oscillations to that contained in radiation now reads: 

osc/ih [^j^j^J VlOOTeV/ V 10.75/ GeY J ^ ^ ' 

The parameter aosci/rh is defined as follows: 



ttosci/rh = min ^1, ■ (32) 

This parameter is most likely 1 if one relies on the constraints on the reheating temperature that result from the 
influence of a moderately massive gravitino on big-bang nucleosynthesis. In models in which the gravitino is very 
massive, ^ 100 TeV, however, such constraints can be evaded. We thus treat aosci/rh as a free parameter. 

Obviously, looking at Eq. ((3T|) . unless Cinf < lO^^Mpi, the modulus is bound to dominate the energy density of the 
Universe at the time of its decay. 

Let us now analyze the constraints presented in Fig. [T] This figure shows the contours of the 6^ry quantity calculated 
numerically, assuming Tjh = 10^ GeV, Hi^i = 10^^ GeV in the left panel and iJinf — 10^ GeV in the right panel, 
— 1. This figure assumes initial (gauge invariant) density contrasts A^'^ = 2 x 10"^ and ^Cinf/cinf = ^finf/(27rcrinf) 
on large scales. The dashed yellow area is excluded by big-bang nucleosynthesis; in order to draw this region, we 
used the results of Ref. Q for a hadronic branching ratio of lO^'^ and initial jet energy 1 TeV. The white dashed line 
indicates the place where n^"* = 0.5 and separates roughly the regions in which either the dark matter or the baryon 
isocurvature mode dominates (see below). There is actually an accidental cancellation of these two modes close to 

that line, it has also been assumed that Co-'^ saturates at 0.5 in order for the numerical computations to proceed 
without errors. In any case, the region in which ^ 1, corresponding to a-mf ^ Hint, is an "unlikely" region, in 
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FIG. 1: Constraints in the {mc,ainf) plane. The yellow shaded region shows the region excluded by the effect of moduli decay 
on big-bang nucleosynthesis. The blue region give the contours of Sm-y for the matter - radiation isocurvature mode (1% and 
10% contours indicated). The white dashed line indicates the place where Sla'' = 0.5. The orange dotted lines indicate the 
standard deviation of the modulus expected from stochastic evolution in its potential during inflation, for new and chaotic 
inflation, as indicated. Left panel: Hint = 10^^ GeV and Trh = lO^GeV. This value corresponds to the energy scale of inflation 
during chaotic inflation. The same value is achieved for small fleld inflation if /i ~ Mpi. In the case of chaotic inflation, 
the standard deviation is given by Eqs. (|A29|I and (|A33|I since the condition < HinfN~^^^ ~ 1.4 x 10^"^GeV is always 
satisfled for rricr < lO^GeV. In the case of small field inflation, the standard deviation is given by Eq. (|A60|I assuming the 
above mentioned condition is also satisfied in this case which amounts requiring that N,j, < lO^'^ . Right panel: _ffint ~ 10^ GeV 
and same reheating temperature. For small fleld inflation, this corresponds to /i ~ 0.22Mpi assuming p = 3. The chaotic 
inflation standard deviation does not appear in this panel because Hinf ~ 10® GeV cannot be realized in this case. For small 
field inflation the standard deviation is again given by Eq. HA60|) which, this time, requires A'^^ < 10''. 



the sense that a-jni is expected to be typically larger than i7inf/(27r) due to quantum effects, see Eqs. ([S]) and ([7]) for 
TOcr,eff < Hinf- The dotted orange lines in Fig.[T]show the standard deviations {Sa'^)^^^ expected from the stochastic 
motion of the inflaton, following Appendix [X] and the formulae Eqs. ([5]) and for two inflationary scenarios and 
two inflationary scales. It is important to stress the following: these standard deviations are measured relatively to 
the "instantaneous" classical value of the modulus field. 

Figure [1] provides a clear example of the power of constraints obtained at the perturbative level on moduli cosmology, 
as the region excluded by the production of isocurvature fluctuations significantly exceeds that constrained by big-bang 
nucleosynthesis. For instance, at a natural scale Too- ~ 1 TeV, the upper bound on crinf (equivalently, on the modulus 
energy density) is more stringent that those from big-bang nucleosynthesis by some two orders of magnitude. This 
figure also clearly shows that taking rria ^ 100 TeV allows to evade the constraints from big-bang nucleosynthesis, but 
not those from cosmological density perturbations if the inflationary scale is large, Hinf ^ 10^^ GeV. One maynote 
that this region is also constrained by the possible overproduction of gravitinos through modulus decay [g^, [H, |69'| . It 
is thus mandatory to require that the modulus potential suffers corrections, in such a way as to reduce considerably the 
modulus energy density at the time of decay, or that the inflationary scale is much lower. Indeed, if iJjnf ^ 10^^ GeV, 
a region devoid of constraints opens up at large values of (Tint and large masses > 100 TeV. The large modulus 
mass then allows to evade the constraints on entropy injection around big-bang nucleosynthesis, while the large value 
of fJinf ensures that the initial modulus - radiation isocurvature fluctuation has become negligible at such a small 
inflationary scale. More specifically, using Eq. pU]) one finds that si'), < (7^ when a > Mpi (i/inf/lO^^ GeV), and 
the isocurvature fraction constrained by microwave background anisotropics is directly proportional to this ratio. 
The constraints obtained in the limit of a high inflationary scale also exclude the possibility of late time entropy 
production through modulus decay. This result is in itself significant as such entropy production is often invoked to 
dilute unwanted relics. 

In order to understand the relative contributions of dark matter - radiation and baryon - radiation isocurvature 
fluctuations, it is useful to break Fig. [1] into two subfigures, each showing one of the two contributions. In the left 
panel of Fig. [2] we plot the same contours as in Fig.[Tl assuming arbitrarily (for the sake of demonstration) C^'' ^ C^'' 
(i.e. setting the baryon isocurvature mode to zero). Similarly, in the right panel of Fig. [51 we plot the contours of (5m7 
assuming ^^"^ ~ (^^•^ (i.e. no dark matter isocurvature mode). These two plots illustrate the complementarity of the 
constraints coming from these two isocurvature modes. It is possible to understand each of them as follows. 
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FIG. 2: Same as the left panel of Fig. [TJ but showing only the contribution of dark matter isocurvature modes in the left panel 
(i.e. setting arbitrarily the baryon isocurvature mode to zero) and the baryon isocurvature mode only in the right panel (setting 
the dark matter isocurvature mode to zero). 



Assuming for the time being, that C^"* ^ C^'' (i.e. no baryon isocurvature mode), the constrained fraction S^-y can 
be written in the following way: 



where 



(34) 



Therefore, the actual quantity that constrains the magnitude of the final dark matter isocurvature mode is x^j. This 
is not surprising since x^-y is, up to a factor Q^/Q,^, exactly equal to the quantity B introduced above. Using Eq. (jl9p . 



one can rewrite x^-y as: 



fin, 1 + V f^x ' + B^n<^ J C« 



It is straightforward to verify that the conditions ([22|) for the generation of isocurvature perturbations lead, if satisfied, 
to a non-zero value of x-^-y. It is however important to note that the magnitude of x^-y also increases with the ratio 

CP/C^'\ One may rewrite the conditions of existence of a dark matter - radiation isocurvature mode as follows, 

neglecting the effect of CctVC-t'"' for clarity: 

VlOOTevJ [lO^Gev) A/pi * VlOOTevJ IvlO^Gevj " ^^^^^ 

Out of simplicity, the dependence on the numbers of degrees of freedom has been omitted in these equations. In 
the above expression, the upper bound comes from the condition fl^"^ <C 1 using Eq. (pij) while the lower bound 
originates from the condition B^fi^"* ^ ^x*^' ^^^^^ formula in Eqs. ((22|) . using the definition of B^^, the fact 

that n<-^ = r2>'T't/T'd (with Tf ~ my^/x{) and Eqs. ((23]) and ([31]). These constraints dehmit a stripe in the 
(cTinf , m^) plane which is consistent with what is observed in Fig. [21 This stripe is actually broader at small values of 
CTinf due to the large ratio C^VCt'' which enhances the modulus perturbations relative to those of dark matter. The 
above formula also neglects the effect of annihilations, which is a good approximation as long as <IC 1, or: 

/^^i"f^//n7r ™- A"'/V V'^V ™X -1/4 



VM] 



PI 



^^"•HlOOTevJ UoSGevi UoOGevJ ^ ' ^'^^^ 
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In practice, annihilations will play a role in suppressing the amount of isocurvature mode in the high mass region 
m„ > 100 TeV and cri„f - 10"^ - lO-^Mpi, as can be checked from the above formulas using Xf ~ 20. 



Turning to the generation of a baryon isocurvature mode, we now assume that Qx 
can be rewritten as: 



(f) 



C,\ \ in which case Eq. ([261 



where 



2xb') 



6 + 



(38) 





'Ai) Ai) 






'd 


) Ai) 
S.7 





(39) 



Therefore the magnitude of the isocurvature mode is proportional to the fraction of energy density of the modulus 
at the time of decay, times the ratio of initial modulus-radiation isocurvature mode to the initial radiation curvature 
perturbation. The baryon constraints thus lie at high values of tiinf , since the production of isocurvature fluctuations 
become dominant when fJ^"* is of order unity. However, as fl^"^ becomes smaller than unity, its weakness can be 

compensated to some level by a large value of Co-'VC-r'^- 



B. Supergravity corrections to the potential 

In supergravity, one expects the potential of the modulus to be lifted by an effective term of the form ztc^H^a^ /2, 
where the factor may change between different eras of the thermal history of the Universe, see Appendix [C] where 
concrete models are discussed. Including such supergravity corrections, the post-inflationary modulus potential may 
be written as follows: 

V{a) ^^mla'±^c'H^<j-<Jof . (40) 

One should recall that we assume the modulus field to be light during inflation, hence the potential ([3D]) refers to 
the post-inflationary epoch only. Moreover, in the following, we denote by "high energy" the regime in which the 
corrections proportional to c^H^ dominate the term proportional to mjcr^ (as mentioned before, although we use 
the expression "high energy", these considerations apply to the post-inflationary epoch only). Then, the quantity c 
controls the mass of the field while ctq represents its minimum, the so-called "high energy minimum" . One may expect 
that the vev at the end of inflation (dinf ) be different from cto since the effective potentials during and after inflation 
a priori differ from one another. We also use the terminology "low energy" to characterize the regime in which the 
corrections c^H^ become sub-dominant and where the potential reduces to rn^a^ /2. Let us recall that we set the 
"low energy minimum" of the potential at u = 0. As a result, the high energy and low energy minima do not coincide 
if (To 7^ (as one should expect on general grounds). 

In order to derive the amplitude of the isocurvature perturbations produced in the supergravity case, we need to 
follow the evolution of the modulus vev and of its perturbations from the end of inflation until the time at which we set 
the initial conditions. In particular, we need to calculate in terms of the inflaton and modulus perturbations during 
inflation. In order to do so, we model the introduction of the effective Hubble mass through a potential that accounts 
for the coupling between the modulus and the inflation, such as V{(j,(f)) = 771^(7^/2 + ni^cfP [l + c^a^ j (3Afp[)] /2, 
which produces the desired +c^H^ effective mass squared. However, the discussion that follows is not restricted to 
this particular potential. In order to study this two field system, we use the formalism of Gordon, Wands, Bassett and 
Maartens [T^I (notice that we have changed notations with respect to Ref. [tO] since, now, a now longer represents 
the adiabatic field but the modulus one). The entropy field is given by the following expression 

Ss = cos 6Sa — sin 95(1) ~ Sa — ^S(j) , (41) 

where we have used that 

cose = ~ 1 , sine = ~ (42) 
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the last equalities following from the test field approximation: a (j). Then, using the perturbed Einstein equation 
for the Bardeen potential on large scales 



&S(T ] , 



one can rewrite the entropy perturbation in terms of Sa and $ only. One obtains 

6s ~ da : — —a^ ^ Sa <i>, 

02/2 



(43) 



(44) 



where we have used that 0^/2 ~ Mpj2 on average. In the above equation, one has neglected a'^ /(p^ in front of 
unity. On large scales, the equation of evolution of the entropy perturbation reads (70j : 



Ss + 3HSs 



(k,s + 39^) 



Ss = 



(45) 



where the expression of Vss can be found in Ref. [70|. In the test field approximation a <^ (j), and assuming powerlaw 
behaviors of a and (p, one can check that Vss — Vaa and 0^ <^ Vac in Eq. so that Ss follows the same quadratic 
equation than cr, and therefore Ssja is conserved. This implies 



Sa 
a 



2a 
'iWa 



S_a_ 
a 



2a 
Zlfa 



inf 



hence, with t = 2/{3H), $ ~ <i>inf and assuming a = initially. 



Sa 
a 



Sa; 



inf 



finf 



ta{t) 



inf 



(46) 



(47) 



This solution happens to match that obtained for a pure time independent modulus potential, see Ref. [4J|. However, 
in the present case, it accounts for the sourcing of the modulus perturbation by the effective Hubble mass. To our 
knowledge, this result has not been obtained before. In the limiting case ^inf 0, one recovers the result obtained 
in Ref. 42] that Sa/a is constant. One can now calculate the modulus - radiation isocurvature fluctuation at the 
initial time, i.e. in the radiation era and after the onset of modulus oscillations, and one recovers the result given in 

Eq. (dni). 

Two additional remarks are in order here. Firstly, in the derivation above, we have neglected the preheating effects. 
This is justified by the following considerations. It turns out that the model investigated here is exactly similar to 



the two field model g 



studied in Ref. 



[71[ with a dimensionless coupling constant given by g — 



Pi)- 



This means that the quantity g^ /{Arn^) = c $ /(24Mp[) <C 1 and that we are never in the "broad resonance" 
regime where preheating effects are important [Tlj . Secondly, we have also numerically integrated the exact equations 
of motion for different cases and have checked that the approximations used above are verified. Above all, we have 
compared the numerical solution for Sa in the post-inflationary epoch to the solution (|47p and have found an almost 
perfect agreement. 



1. Case +c^H'^/2 

As mentioned previously, the details of the calculations that follow can be found in Appendix |B1 In order to 
re-compute ilo-^osci, one must follow the evolution of a between the end of inflation (where a = ai^f) and H — mu in 
the situation where the potential is dominated by the Hubble scale corrections. This evolution is characterized by c, 
(To and p, the latter defining the evolution of the scale factor: a cx fP , so that H = p/t. These parameters enter in the 
following combinations: 

^^3p_3^ = (/i+l)^-pV, (48) 

Depending on the magnitude of c, v may be real or imaginary, which gives rise to different evolutions. We examine 
each of these cases in turn. 
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Let us first start with the case c < (p + l)/p (real ly). In Appendix [Bl it is shown that f^CT^osci can be expressed as 
follows, see Eqs. (|BT9)) and ((B20)) 

1 1 /_ \ 2(/j.MD + l — ^^a-id) 

^^.o.. ^ r ^) + i^? f ^-^) r ^) , T.. > T,. , (49) 



6 ' \MpiJ 6 ' \ Mpi J \ Hini 

The quantities indexed with "r,d" (resp. "md") refer to the radiation dominated (resp. matter dominated) era. These 
equations are the supergravity counterparts of Eq. ([50)1 in the case where v is real. In the above expressions, the first 
term is the contribution originating from the particular solution of the equation of motion while the second term is 
due to the homogeneous solution. A crucial difference between the scalings of the homogeneous and the particular 
solution is the redshift factor (m„/7Ji„f)2(A'MD+i-''MD) ^ i foj. ^j^g latter. The prefactors Ai, A2, Bi and B2 are all 
of order unity, see Appendix |B1 The quantity ria-,osci will be dominated by the homogeneous solution contribution 
whenever: 

ki„f-ao|(^|^j » kol , (51) 

assuming for simplicity Tih < Tosci which is the most generic situation. 

Let us now discuss the cosmological consequences for the two cases in which the particular and the homogeneous 
solution dominates at late times, starting with the former case. Then, the constraints in modulus parameter space are 
straightforward to derive. Using the results obtained for the purely quadratic case, for which f2cr,osci = (o'inf/Mpi)^/6, 
one can put an upper limit on cto since, in the case considered here, r2o-,osci has the same form, crinf being simply 
replaced with ctq . In Eq. (j36[) . we have established the conditions for the existence of a dark-matter isocurvature 
mode. It is clear that if ao [fiinf in Eq. ()36|1 ] is smaller than the lower bound, then there is no isocurvature mode 
and the scenario is compatible with the cosmic microwave background data. For X[ ~ 20, TOo- = 10^ GeV and 
= 10^ GeV, this gives 

ao < 10"^"Mpi . (52) 

From previous analytical calculations, one expect this bound to scale with the parameters as H^^^T'^^^'^mV^ , see 
Eq. (|36p replacing (Tinf by ctq. However, such a value remains well below the typical displacement expected from the 
quantum jumps of the modulus in its potential during inflation, unless iJinf ^ 10^ GeV (see Eq. [5]). 

This brings us to the other extreme case, in which the homogeneous solutions dominates the evolution at late times. 
As a clear example of this situation, consider ctq = 0, but Uinf 0. Physically, this corresponds to the situation 
in which the minima of the effective potential after inflation coincides at high {Hini > H > rria) and low energy 
{rria > H). The vev Cinf is here non-zero, either because the effective minimum during inflation does not coincide 
with that at latter times, or because Uinf is subject to quantum effects. We thus treat Cinf as a free parameter as 
before, and one obtains the constraints in the modulus parameter space presented in Fig. [3] 

This figure clearly confirms the power of constraints obtained at the perturbative level. With respect to Fig. [1] 
which corresponds to the purely quadratic case, one can see that the constraints from big-bang nucleosynthesis and 
from the production of isocurvature perturbations have moved towards higher values of Cinf . This is understood easily: 
in the purely quadratic case, the modulus energy density remains constant between the end of infiation and the onset 
of oscillations, while it decreases in the present case. The white dashed line, which indicates the locus of 0^=' = 0.5, 
serves to delimit the constraints derived from dark matter - radiation and from baryon - radiation isocurvature modes. 

Nevertheless, the constraints obtained still preclude the possibility of having a heavy modulus with an arbitrary vev 
at the end of inflation. In particular, even if the effective minima of the potential during inflation coincides with that 
lower energy, jcTint — ctqI should depart from zero by the value expected from quantum motion during inflation. In the 
left panel, for -ffi„f — 10^^ GeV, it is found that the constraints from the production of isocurvature fluctuations extend 
significantly below this bound, hence there is no apparent solution to the moduli problem. As i?inf becomes much 
smaller, some region of parameter space opens up at large values of dinf and large values of m„ as in the quadratic case, 
and the typical stochastic displacement also decreases. For iJjnf — 10^ GeV, corresponding to the right panel of Fig. [31 
there is however little room between this lower limit for <Ti„f and the region excluded by isocurvature fluctuations. 

It is noteworthy to recall that the situation depicted in the above figure is realized by two concrete models of 
infiation discussed in Appendix [Cl 

Let us now turn to the case c > (/i + l)/p (imaginary v). For p = 2/3, this case corresponds to c > 3/4. Writing 
u = iv, one has v > Q growing with c. In this case the energy density stored in the modulus at the onset of oscillations 
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FIG. 3: Same as Fig. [T] but considering positive squared Hubble scale corrections to the modulus potential, with (? — 0.5 and 
CTo ~ 0. The parameter v is given by v^j^ ~ 0.027 and ~ —0.063. For the reheating temperature considered in this figure, 
one always has Tosci > Tj-h and, hence, only the value of i^md matters, see Eq. (|49[). Therefore, one is indeed in the case v > 
Left panel: Hi^i = lO" GeV, T^h = 10^ GeV. Right panel: ffi„f = lO'' GeV, = 10^ GeV. This case is only relevant for small 
field inflation. The standard deviations due to stochastic motion are calculated as in Fig. [1] 



can be expressed as, see Appendix [Bl especially Eqs. (|B34I) and (|B35p 
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(53) 
(54) 



These equations are the counterparts of Eqs. (|49)) and ((50)) in the case where the quantity v is complex. The coefficients 
^3 and Ai are defined in Eq. (jB36|) . The prcfactors and are defined in Eqs. (|B37|) and (IB38|) and are of order 
unity. The difference with the case v > Q (with prefactors Ai and A2) comes from the fact that, now, the numerical 
prefactors A^ and Aa may become quite small if v is pure imaginary and its modulus is large. The prefactor ^3 
indeed scales as: 



A. 



27r 



Maid +2 



-i>MD7r/2— mmd — 3 



^mdI > 1 ■ 



(55) 



Hence, as c grows beyond 1, the value of ricr,part,osci decreases exponentially. As argued by Linde |73 |. this could 
alleviate the moduli problem, although it is notoriously difficult to construct explicit models in which > 10. 

Let us now study in more details the physical consequences of the above expressions. Ignoring the factors of 
order one for simplicity, and using the explicit expression of the coefficient A^, see Eq. (jB36p . the particular solution 
dominates whenever the following condition is valid 



Finf — CTqI 



« koir 



V 2 



(56) 



assuming for simplicity Tosci > ?rh. Then, one finds that the bound on do obtained previously in the case of a real v 
for i?int ~ 10^'^ GeV is now loosened by: 



-10 



2(mmd+2) 



-i>MD7r — 2(mmd+3) 



PI 



(57) 



for ma — 10^ GeV, scaling approximately as raj'^ . To provide concrete estimates, for (? = 10 (j^ ~ 2.04 assuming 
matter domination), the r.h.s. becomes 2 x 10~^Afpi; for (? = 20 (i^ ~ 2.94), it is 1.2 x 10~®Mpi. Interestingly, 
these values always remain smaller or are at most (for (? = 20) comparable to the standard deviation expected 
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FIG. 4: Same as Fig.[T] but for a modulus potential receiving a positive Hubble mass squared correction with c? = 2, and ao — 0. 
This implies ul^^ ~ -0.64 and ~ -0.44. Left panel: Jfinf = lO" GeV, T^h = lO'' GeV. Right panel: = lO'' GeV, 

Trh = 10^ GeV. For this value of the reheating temperature, one always have Tosci > Trh in the parameter space considered in 
this figure. 



from stochastic motion of the modulus if i?inf ~ lO'^'^ GeV. Therefore, solving the moduli problem in this way would 
require not only a large value of in order to lessen the modulus energy density, hence the transfer of isocurvature 
perturbations, but also a small inflationary scale <C 10^^ GeV in order to diminish the magnitude of stochastic motion. 
Note also that, for smaller values of -ffinf, some parameter space opens up at large values of Cinf and large mg-, as 
discussed before. 

If, on the contrary, the homogeneous solution dominates, then the situation is slightly different and the constraints 
in the modulus parameter space are presented in Fig. 3) The comparison with Fig. [31 which presented the constraints 
for real v (and ctq = 0) is straightforward. One still obtains an upper bound on jfjinf — ctqI if -ffinf ^ 10^^ GeV, albeit 
displaced to larger values due to the fact that Q.^, now scales as H^^^^ instead of //m+i-2i>^ j ^i^^ energy density 
contained in the modulus decreases faster. This upper bound can be written as: 

kinf-fTol < 3x lO-^Mpi , (58) 

with the following set of parameters: i?inf = 10^'^ GeV, Trh = 10^ GeV, = 2 (z/ ~ 0.8 assuming matter domination 
after inflation), and the prefactor has been evaluated for nia = 100 TeV. As before, one flnds that this bound scales 
as HlJ^^mP^T~^^^ . However, even for = 2, this upper bound is smaller than the amplitude of stochastic motion of 
the modulus in its potential during inflation, so that this cannot be considered as a viable solution. One therefore 
has to require i/i„f ^ 10^'^ GeV, in which case some parameter space opens up at both small values of a-^i (because 
the amplitude of stochastic motion is smaller at smaller 7?inf) and at large values of CTinf (where the isocurvature 
fluctuation becomes much smaller than the curvature perturbation). 

Finally, let us end this section with the following remark. Above we have discussed the cases of v real or pure 
imaginary in both the radiation dominated and the matter dominated epochs. Of course, there are also two mixed 
cases corresponding to a real u during the matter dominated era and a pure imaginary v in the radiation dominated 
era and the opposite situation. This would be relevant for the discussion above only if Tj-h > Tosci, since the matchings 
would have to be modified. Since this is not the most generic case, we do not consider this situation in this paper. 



2. Case 

If the modulus receives a negative Hubble mass squared contribution after inflation, the minimum of its potential is 
destabilized and as a result, the modulus will move until this negative contribution is balanced by higher order terms 
in the potential, cx or even non-renormalizable contributions. Let us assume for instance that the next term in the 
modulus potential takes the form: 
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Then, as shown in Ref. [42] the effective potential may be approximated at high energy H ^ ma/c by: 

T//' N 1~2tt2/ \2 , -^ri (ct — Cr„) 

Via) ^-cH [o^ .„) + (^^,), , (60) 



where the time-dependent quantity (T„ can be expressed as 



(n + 3) 



I -I l/(n+2) 



iJ^M^i 

An 



(61) 



and = (n + 2)c^. The value (t„ corresponds to the local minimum of the potential. Constant terms in V{a) 
have been omitted in Eq. (|60p . as well as subleading terms when compared to the last term on the right hand side, 
see Ref. [13]. One crucial difference between this effective potential and that obtained for +c^H^ is the fact that the 
local minimum now depends on H and thus evolves in time. 

One should distinguish three phases of evolution depending on which term in the potential dominates. If a-mf 3> 
cr„jnf = crn{H = i?inf), then the field initially evolves in the high order part of the potential given by V{<t) ~ 
A„cr^+"/[('^ + 4)!Mp[]. As shown in Appendix iBl the field is then driven to cr„ within a fraction of e-fold of order 
(finf/o'n^inf)^^"'''^-'^^. At this stagc, the effective c^H^{a — of tl^c potential takes over the high order part. In 
order to model this case (Tinf 3> cr„^inf , we simply assume that, starting from the end of inflation, the field evolves in the 
{a — an)^ part with an initial value of order i7n,inf- In this situation, we also neglect the initial kinetic energy since 
kinetic energy is strongly damped when the field evolves in the high order part of the potential. Furthermore, this 
is conservative in the sense that it underestimates the energy density contained in modulus oscillations at late times, 
and therefore tends to loosen slightly the constraints derived. If (Tinf ^ Cn^nf, then one should directly approximate 
the potential by (?H'^{a — ct„)^ with Cinf as the initial condition. All in all, it suffices to solve for the evolution of a 
in the potential c^H'^{a — cr„)^ with an initial condition (Teff.min = niin(crinf , cr„^inf). 

Thus ignoring the high order part of the potential, at high energy H ^ ma/c, the solution for a comprises a 
solution to the homogeneous equation crhom and a particular solution (Tpart, as before, see Appendix iBl The scaling 
of the homogeneous solution is similar to that obtained in the previous section with a +c'^H'^ effective squared mass 
term, but the particular solution scales differently due to the time dependence of a„ . As shown in Appendix [B] 
Cpart o(. cr„ with a prefactor q;„ of order unity. The explicit expression of a„ can be found in Appendix [B] see 
Eq. (|B43p . Then Cpart c>c t^2/(n+2) t}igi± the energy density of the modulus associated to this particular solution 
scales as Per, part H^a^^j.^ (x ij2("+4)/(n+2)^ Consequently, the particular solution contribution to the energy density 
at the onset of oscillations can be written as 



where we have ignored all factors of order one. In the above expression, (Tpart, inf represents the initial value of the 
particular solution at the beginning of the era driven by c^H^{(t — (t„)-^, that is to say 



^part,inf — O^n^^n, inf — 



jn + Sy. 2 2 ..n 
T C i?infAipi 



l/(n+2) 

<//"+2)a<("+2) . (63) 



Concerning the homogeneous solution, its evolution is the same as in the previous section (see also Appendix |B]) . 
This means that one should again distinguish the case where v is real or imaginary and should treat separately the 
situation where the onset of oscillations occurs before or after the reheating. Straightforward calculations, similar to 
the ones already performed in the previous sections, lead to our final expression of r^o-.osci 



1 / \ 2 / \ 4/(n+2) -I / \ 2 / \ 2[;iMD+l-Ro(i^MD)] 

1 / tTpart,inf \ / »V \ \ ( ~ Crpart,inf \ / ^ ^'^ 

1 / fpart 



> T,h , (64) 



6 I A/pi J V H, 



2 / X 4/(n+2) 



1 /_ _ \2 / ™ \ 2[^iRD + l-Ro(mD)] / TT 

J_ / Crcff,inf — Crpart,inf \ / "T-ct \ / J^rh 

6 ^ Mp^ J [h^J \H~i 



2[/tMD+l-Ro(l/MD)] 



< r,h , (65) 



where we have neglected the factors of order one and where the appearance of the real part of v accounts for both 
possibilities real or imaginary). 
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FIG. 5: Same as Fig. \T\ for a potential receiving a negative Hubble mass squared correction with = 0.5, and n — 0. 
Hint — 10^^ GeV (for Hinf = 10® GeV, constraints from isocurvature fluctuations and from big-bang nucleosynthesis vanish). 

When discussing the case +c^H'^, we considered two cases, one in which the particular solution dominates, the 
other in which the homogeneous solutions dominates. We cannot do so here, because the particular solution is 
entirely determined by Hmi and no longer dependent on the magnitude of cto. Both contributions have to be considered 
together. Furthermore, we recall that (Toff^inf = min(crinf , cr„jnf )■ Therefore, whether the particular or the homogeneous 
solution dominates in Eqs. ([M)) and (p5|) depends on n, /z, u and m^/Hiai (assuming Tosci > T^h). 

Regarding the fluctuations of a, it is not possible to follow analytically the evolution of Sa/a from the end of 
inflation until H ~ m„ due to the non-linearity of the potential. For the sake of the argument, we thus assume that 
the initial conditions arc the same as in the previous cases, namely Eqs. (Pjl. pH]) . 

Let us flrst discuss the case n = and assume for the sake of discussion that — 0.5. Then, j/md is imaginary since 
(? = I. Furthermore, 2(^md + 1) = 1 and 4/(n-|-2) = 2 in this case, so that the particular solution is always negligible 
in front of the homogeneous solution. The contribution of this latter to the energy density is nevertheless suppressed 
by rricr/Hi^f and fl'^'^ is so small that the constraints from big-bang nucleosynthesis are signiflcantly weakened, see 
Fig. El they now allow moduli masses above SOOGeV for all tiinf. This constraint does not depend on crjnf because 
ri<<i hardly depends on crinf since (f7off,inf - crpnit.inif ~ o-part,inf both hmits dint < crpart,inf and ai^i > o-part,inf • 
However, the constraints from the production of isocurvature fluctuations are quite signiflcant in this case because 
the modulus can perturb signiflcantly the perturbation spectrum of dark matter eventhough it does not dominate the 
energy density at the time of decay. Strictly speaking, this is true as long as Hinf ~ 10^^ GeV, as these isocurvature 
constraints are less stringent for Hint — 10^ GeV. 

When n > 1, r^o-^osci becomes large enough to produce substantial isocurvature fluctuations in nearly all of the 
parameter space, see Fig. [51 For n = 1, the homogeneous solution still dominates over the particular solution. 
The constraints can be derived from the discussion of Sec. 1111 B II (i.e. +(?H^ potential, in the case in which the 
homogeneous solution dominates) provided one replaces (dinf — ctq)^ with (cinf — cTpart.inf)^ — fpart inf ' which is 
a function of n and ffinf, but which does not depend on dint- Accordingly, $1^'^ is sufficiently large to produce 
signiflcant isocurvature perturbations because of the scaling of CTpart.inf with n: for n = 0, CTpart.inf ~ ^^inf, but for 
n > 1, CTpart.inf ^ ^^inf, scc Eq. ((63)) . As before, we find that the isocurvature mode becomes small enough to satisfy 
the constraints from cosmic microwave background data if Hi^t ^ 10^"^ GeV and a\nf ^ Mp\. 

C. Modulus production by inflaton decay 

The above discussion has implicitly assumed that no modulus was produced after inflation. However it seems 
reasonable to assume that the inflaton can decay into the modulus sector, through possibly Planck suppressed inter- 
actions. We take a branching ratio such that each inflaton produces N„ moduli. Unless the modulus and the inflaton 
are coupled one to the other, one should expect A^^ ^ l/5*,rh, which means that at most, moduli are produced at the 
same rate than other light particles. If the inflaton is more strongly coupled to the visible sector than to the modulus 
sector, one should expect a much lower value of N^,. 

Since the modulus mass is generically much smaller than the inflaton mass, the moduli produced through inflaton 
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decay are ultra-relativistic, with energy E^^^ ~ These particles do not thermahze but redshift to non- 

relativistic velocities. As their momentum redshifts away according to Pa = m^/2(ai.h/a), the particles become non 
relativistic, when p or Oih/an-rci = '2.m„/m^. The corresponding temperature T,i_roi is given by 

1/3 1/3^1 2TOo- 

5*„-rcl^n-rcl = g„ijr\.—— ■ (66) 

For temperatures T > T^-rci, Pa^4> oc 1/a^ i.e. it scales as radiation while for T < Tn-ioi, it scales as a pressureless 
fluid, Pa-^4> oc Just after reheating, the energy density of the moduli that were produced by inflaton decay is 

given by 

rh AT rh /'C'7^ 

If the temperature Tn-roi is smaller than Td, then p^^^ will scale as the radiation until the modulus decay. As a 
consequence 



Pi 



<d 2 



r„-,ei < Td . (68) 



If, on the contrary, Tn_i-ei > 7d, the energy density Pa^4> will then increase with respect to that of radiation in the 
era following T^-j-d and preceding modulus decay. Therefore, the ratio of inflaton produced moduli energy density to 
radiation energy density immediately before modulus decay can be written as: 



p^^ 2 Td 

- "^-^HlO-V U*,n-rcJ V10.75J VlOOTevJ VlOi3GevJ GeV j ' ' > 



(69) 

where one has used the expression ([3|) of Tj. Let us notice that, in most of parameter space, Tn-roi > 2d- 

The next step is to study how the previous considerations will impact the calculations developed in the previous 
sections. Obviously the amount of energy density stored in the modulus oscillations remains unchanged. However the 
amount of isocurvature fluctuation produced at a same value of fJ^'' is reduced 9]. Indeed, inflaton produced moduli 
inherit the same spectrum of perturbations than radiation, and therefore there is no initial isocurvature perturbation 
between those moduli and the radiation fluid. Including these inflaton produced moduli, the initial value of the 
modulus curvature perturbation now reads: 

C«' =7C« + (l-7)d^ (70) 
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The quantities in this equation should be understood as follows: and correspond to the curvature perturbations 

of the modulus and the inflaton acquired during inflation, as before. The curvature perturbation (^i'^ ' is the curvature 
perturbation for the modulus that should be used in Eqs. ([55)) and ([55)1 for the calculation of the final isocurvature 
perturbations. Finally, 7 denotes the ratio of the energy density of moduli initially present at decay to the total 
amount of moduli at decay (those initially present together with the inflaton produced moduli). The final effect 
is to modify the initial modulus - radiation isocurvature perturbation by a factor 7 in Eqs. ()35p and ()39p . namely 
— > 72;^^. It is also useful to define 7' = p'^^fj,/ P^"^ (with p^"^ the amount of energy density stored in the modulus 
oscillations immediately before decay): 

(71) 



7' + l 



so that, if Ta-i-ci > T^: 

^' " iw^) n7-. (lo^) (lo^) ^ ' > ■ 



1/3 



In order to obtain this formula, we have used Eq. ([M)) in order to express p^^^ and Eq. ()25p to express p^<^ [or, 
alternatively, Eq. ()3ip with the term cr?„f/Mp[ replaced by Orio-^osci in order not to be restricted to the quadratic case]. 
Notice that we have taken ttosci/rh = 1- 

If, on the contrary, Tn-jci < 7d, then the 7' factor can be expressed as 

y.oa..«-"(j^)^(|^)-(^)-(^)-', ,7., 

where, this time, p^^^ has been obtained from Eq. (|68p . In these equations, r^o-^osci should be understood as corre- 
sponding to the oscillations of the modulus. It docs not include, in particular, the moduli produced through inflaton 
decay. 

To study the effect of such modulus production through inflaton decay, we first assume that the potential is 
purely quadratic, i.e. c = 0. The results are presented in Fig. [Tj Compared to Fig. [1] one finds that the big-bang 
nucleosynthesis constraints now exclude all moduli masses below 100 TeV. This is expected insofar as the amount of 
moduli energy density produced through inflaton decay is sufficient to disrupt big-bang nucleosynthesis; since this 
amount does not depend on CTinf , contrary to the amount of energy density stored in moduli oscillations, the big-bang 
nucleosynthesis constraints also do not depend on aiai- Note that the inflaton may also decay into gravitinos, with 
similar consequences for big-bang nucleosynthesis, see [73j . 

The constraints from isocurvature fluctuations are pushed to larger values of CTinf, since the factor 7 becomes small 
when the energy density produced through inflaton decay far exceeds that stored in modulus oscillations. Conversely, 
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FIG. 8: Same as Fig. |4l for a receiving a positive squared mass Hubble contribution with c? — 2, but now accounting for 
inflaton produced moduli with Na — 10~^ and an inflaton mass = 10^'^ GeV. Left panel: Hi^f = 10^^ GeV, Trh ~ lO'' GeV. 
Right panel: Hi^t = lO" GeV, T^h = 10^ GeV. The stochastic effects are calculated with C; —2. 



a larger value of (Tinf yields a larger value of rio-,osci hence a larger value of 7. Assuming for simplicity Tn-rci > 
one can check that the suppression of isocurvature fluctuations becomes effective for 

In order to obtain this expression, we have used Eq. (I72p and have written the condition 7' 3> 1 (which is equivalent 
to 7 <C 1) in the quadratic case, namely llosci ^ '^tnf I ^'-^fv This allows to understand, at least qualitatively, the trend 
shown in Fig. [71 

The effect in the case where the potential receives supergravity corrections is rather straightforward to guess: big- 
bang nucleosynthesis constraints remain unchanged as compared to the above Fig. [71 but the contours depicting the 
amount of isocurvature fluctuations produced are shifted toward higher values of result of the redshifting 

of the energy density stored in modulus oscillations after inflation, yielding a smaller value for 7. For instance, 
considering the case = 2 as in Fig. [H one obtains the constraints depicted in Fig. [51 

Quite interestingly, the production of moduli through inflaton decay, while reducing the overall amount of isocurva- 
ture fluctuation, does not allow to find a solution to the moduli problem with arbitrarily high [~ ©(Mpi) in particular] 
vev at the end of inflation if ffinf — 10^'^ GeV. Strictly speaking, the amount of isocurvature fluctuations produced 
for this value of the Hubble constant and crinf ~ 0.01 — 1 M-p\ is of the order of a few percent, therefore it is not 
excluded by present cosmic microwave background data. If present at this level, it could actually be detected by 
the upcoming generation of instruments. One should also emphasize that we have considered a rather conservative 
case, in the sense that the modulus is produced at a comparable rate than other light particles; if is decreased, 
the amount of isocurvature fluctuations would increase. At smaller values of i?inf , the isocurvature constraints have 
disappeared, due to the combined effect of the partial erasure associated with moduli production in inflaton decay and 
a smaller initial isocurvature fluctuation. Finally, independently of i?inf, the production of moduli through inflaton 
decay significantly worsens the effect of moduli on big-bang nucleosynthesis. The success of big-bang nucleosynthesis 
now requires both mo- > 100 TeV for all CTinf- 



IV. SUMMARY AND CONCLUSIONS 



Let us first summarize the results obtained. We have shown that the decay of a generic modulus tends to produce 
strong isocurvature fluctuations between dark matter and radiation, or between baryons and radiation. The amount of 
isocurvature fluctuations produced, relatively to the total curvature perturbation, depends on several parameters: the 
value of the initial modulus - radiation isocurvature perturbation (in units of the total initial curvature perturbation), 
and the amount of energy density stored in the modulus oscillations at the time of decay, relatively to that contained 
in radiation, in particular. We have discussed in some detail the evolution of the modulus energy density and of its 
perturbations from the end of inflation onwards, for a variety of possible moduli effective potentials, assuming the 
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modulus remains light during inflation. We have then translated the constraints derived from the analysis of cosmic 
microwave background data into constraints in the modulus parameter space mo- — CTinf • We find that the constraints 
associated with the production of isocurvature fluctuations significantly exceed those from big-bang nucleosynthesis 
in this parameter space. One reason why the constraints obtained cover most of the rria — crinf parameter space is 
that the modulus will produce a large baryon isocurvature mode if ri^"* ~ 1, while it will produce a dark matter - 
radiation isocurvature fluctuation if il^'' <C 1 (but ri^^* large enough to affect the dark matter). 

Evading the constraints from big-bang nucleosynthesis and from the generation of isocurvature fluctuations requires 
one of the following conditions to be satisfied. 

Firstly, if the modulus potential is time independent (i.e. it does not receive supergravity contributions at any time), 
it is mandatory that: (i) -ffinf <SC 10^ GeV and the modulus initially lies very close to the minimum of its potential, 
within ~ 10~^°Mpi (depending on the modulus mass, here taken to be ~ ITeV, see Fig. [T]); or (ii) m„ > 100 TeV, 
<^inf ^ (^^inf/lO^'^ GeV) Mpi and Hinf ^ 10^'^GeV. Solution (i) ensures that the moduh energy density at the time 
of its decay is sufficiently small to affect neither big-bang nucleosynthesis nor the dark matter perturbations (which 
turn out to be more sensitive probes than the former in this region of parameter space) . The constraint on the Hubble 
parameter corresponds to the requirement that the stochastic motion of the modulus during inflation remains small 
enough as compared to the bound on the final effective displacement of the modulus vev. Such a solution might be 
realized in models in which the modulus is bound to remain close to an enhanced symmetry point in moduli space due 
to the friction caused by its coupling to other light degrees of freedom [Ijlzl, (zH]. Solution (ii) is typical of particle 
physics models which achieve a high modulus mass scale. The constraints on a-mi and most particularly the bound 
on ifinf directly come from the constraints due to isocurvature fluctuations. If Hint <C 10^'^ GeV as required by this 
solution, the tensor modes should be unobservable by upcoming cosmic microwave background missions. 

Secondly, if one assumes that the modulus potential receives supergravity corrections after inflation (but remains 
light during inflation, as realized in some models of inflation discussed in Appendix [Cjl . one has to require Hmf <C 
10^^ GeV if the minima of the effective modulus potential at high energy {H 3> nicr) matches that at low energy 
{H <C nicr)- Then the moduli problem can be solved either for large tJinf ~ Mpi, provided > 100 TeV, or for 
vanishingly small dinf (in units of Mpi, see the main text for details) and arbitrary TOo.. If the minima do not coincide, 
the situation is very similar to that discussed for a time independent modulus potential. 

Finally, we have recalled that if the modulus is heavy during inflation, cither because i?inf ^ rn^, or because the 
modulus receives an effective Hubble mass during inflation, then the isocurvature fluctuation between the inflaton and 
the modulus disappears (Blj . This provides a natural solution for the moduli problem at the perturbative level, but 
it does not automatically satisfies the constraints on entropy injection at the time of big-bang nucleosynthesis. First 
of all, either ma- ^ 100 TeV, or the minima of the effective modulus potentials during inflation, after inflation and at 
low energy (H <^ rria) match one another. In this latter case, one also needs to ensure that the quantum jumps of the 
modulus in its potential does not result in too large an energy density at the time of big-bang nucleosynthesis. This can 
be done using the discussion of Appendix El which discusses the stochastic motion of the modulus in the presence of a 
Hubble effective mass CiH. Strictly speaking the results are valid in the slow-roll regime, which requires Ci < 3/4, but 
they should remain valid as long as as Cj is of order unity. For chaotic inflation, the standard deviation {Sa'^)^^^ is given 

by the Bunch-Davies expression (|A40p . This leads to the bound of Eq. (|A44p which gives ((5a^) ~ 3.0 x 10~^Afpi 
(here for cf = 0.5). For small field inflation, the Bunch-Davies expression is given by Eq. (jA74p . In the case where 

/i ~ Mp\, which is the case if Hini — 10^^ GeV, then one can also use Eq. (|A77p and this leads to ((5cr^) — lO^^Mpi. 
These values are in conflict with the big-bang nucleosynthesis constraints if ^ 1 TeV and the modulus potential 
becomes time independent after inflation, since these constraints require a-mi < 10^^" Mpi (see Fig.[T]). When turned 
into a bound on the Hubble constant during inflation, this becomes Hini ^ 10^ GeV, a low value indeed. Even if 
the modulus potential receives supergravity corrections after inflation, there is a potential conflict, since the big-bang 
nucleosynthesis constraints impose Cinf < 10~^ Mpi for = 0.5, but Uinf < 10~^ Mpi for = 2, see Figs. [31 SI 

We have also considered the consequences of a possible production of moduli through inflaton decay. Since those 
moduli inherit the fluctuations of the inflaton, which correspond to those of radiation, this tends to reduce the initial 
modulus - radiation isocurvature fluctuation, all things being equal. However, it also aggravates the effect of moduli 
on big-bang nucleosynthesis, and one generically finds that all moduli must have masses nia- ^ 100 TeV independently 
of cTjiif . Furthermore, one must still require that the effective minima of the modulus potential during inflation, after 
inflation and at low energy match one another, or Hinf ^ lO^'^ GeV as above. 

In short, the moduli problem is worse at the perturbative level. A clear trend emerges from the above calculation, 
namely the lower the inflationary scale, the easier it is to solve the moduli problem. A significant amount of late time 
entropy production could alleviate the moduli problem, provided baryon isocurvature perturbations are not produced 
at the decay of the entropy producing component. As mentioned above, one needs to have a reheating temperature 
(after entropy production) higher than the baryogenesis scale; alternatively the fluctuations carried by the entropy 
producing fluid could be similar to those carried by the baryons. Those are certainly non trivial constraints. 
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APPENDIX A: EXPECTED VALUE FOR a-.^t 

We now analyze the stochastic (quantum) behavior of the modulus field in the case where the total energy density is 
still dominated by the vacuum energy of the inflaton field. A similar analysis was performed in Refs. [sl but focusing on 
the chaotic inflationary scenario and for negligible initial values of the modulus field. Moreover, the stochastic nature 
of the inflaton field was also ignored. Here, we relax these assumptions and generalize the results of Refs. [ol. [5l|. 

The problem treated here bears close resemblance with the problem tackled in Ref. fl^ where the behavior of the 
quantum quintessence field during inflation was studied (the quantum behavior of the infiaton field being also taken 
into account). For this reason, we will follow a similar treatment. 

According to the formalism of stochastic inflation, the coarse-grained inflaton field (j> obeys the following Langevin 
equation [H, |4§| 

where is a white-noise field such that {^4>{t)S,^{t')) = d{t — t') and where a prime denotes a derivative with respect 
to the field. The stochastic evolution of the modulus field a is also controlled by a Langevin equation which, in the 
slow-roll approximation, reads 

where ^o- is another white-noise field such that 

{Ut)Ut'))^5{t^t'), {Um4t'))^0. (A3) 



The solution of the Langevin equation (|A2p depends explicitly on but also on the inflaton noise through the 
coarse-grained field (/>. Since the modulus is considered as a test field, H only depends on cf), hence all primes superscript 
mean differentiation with respect to (p; obviously, all primes superscript on (resp. Va) denote differentiation with 
respect to (j> (resp. a). 

In order to find an approximate solution to Eqs. (|A1[) and (jA2p . one may try to use the same perturbative technique 
as the one used in Refs. [t^, llHl- Therefore, we expand the inflaton and modulus fields about their classical solution 
and write 

4>it) = c/jciit) + S(j)i + 6(1)2 + ■■■ , a{t) = adit) + Sai + 5(72 + ■■■ , (A4) 

where the first terms in the expansions are linear in the noise, the second ones are quadratic in the noise and so on. 
As mentioned before, the approximation made in Ref. Q consists in ignoring the stochastic nature of the infiaton 
field, (f> = (f)ci- 

Let us first examine the solution for the inflaton fleld. Expandin g up to second order in the equation of motion, we 
obtain two linear differential equations for d(j)i and S(j)2, see Refs. [771. l78l|. namely 

and 

^ + 2Ml,H"{(lyd)H2 = -Af|ii?"'(0ci)^0? + |^i?i/2((/.ci)i/' (001)^01^0(0 • (A6) 

These equations can be solved by varying the integration constant. Let us flrst consider the equation for 5(j)i. If the 
initial conditions are such that S(f>i(t = ti^i) = 0, then the solution reads 
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We are now in a position where the various correlation functions can be calculated exactly. Since S(j)i is linear in the 
noise ^, the mean value obviously vanishes {S(f>i) — 0. The two-point correlation function for S(f>i can be calculated 
as: 

{SMti)SMt2)) = -^H'{h)H'{h) dr-^^ . (A8) 

Once the inflaton potential has been specified, it is more convenient to carry out these integrals using the classical 
value 4>c\ as a variable instead of t, thanks to the relation derived from the classical trajectory: 

d0ci = -2H'Mlydt . (A9) 

For instance, the two-point correlation function calculated at the same time, i.e. the variance, reads [ttI . pTSj 

^jy ■ 

Detailed calculations of these integrals will be given below for various prototypical models of inflation. 

We now turn to the equation of motion for the second order perturbation 54)2 ■ It can be solved by following exactly 
the steps that were described before. Then, the solution can be written as [73, pTSj 



-H'Ml,,j^ dr—Scflir) + ^ drH'/'SMr)Ui^) ■ (All) 



As expected the second order perturbation is quadratic in the noise. One can easily evaluate the mean value of 
(5(/)2(i), taking into account a factor 1/2 which originates from the fact that the Dirac ^- function appearing in the 
noise correlation function is centered on an integration limit, see Refs. [ttI pTSj 

Let us now turn to the modulus case when it has a sufficiently flat potential to acquire an independent quantum 
noise besides the inflaton. As for the inflaton case, it is easy to establish that the equations of motion for the perturbed 
quantities Sai and Sa2 are given by the following expressions [t^ 

d5a2 ^ C(^ci) , Fj(aci)g'(0ci) ,, ^ K(^ci)g"(0ci) ,,2 V;(c7ci)g^^(0ci) , ,2 

+ 3g^(0el) '^'^^'^"^ " ^^^"^ + '^'^^^^ • ^^'^^ 

Although these equations look quite complicated, they can be solved easily because (by definition) they are linear. 
Assuming that the modulus potential does not depend explicitly on time (that is, other by its dependence on cr), the 
solution for Sai reads: 



— ^— — d^l r -I- — — riaij) 



dr , (A14) 



and, as required, is linear both in the quintessence noise and (through 5<j}i) in the inflaton noise The above 
formula is not valid anymore if Va contains an explicit dependence on time, for instance if the modulus mass receives 
Hubble term corrections. In this particular case, which will be discussed further below, one has to extract the explicit 
time dependence out of the potential, then proceed as above. Unless otherwise said, we assume in the following that 
Va does not contain any such explicit dependence on time. 

As is obvious, 5ai has a vanishing mean value, (i5cri) = 0, but a non-vanishing variance given by the sum of two 
contributions originating from the inflaton and quintessence noise variances, namely j76| 

(A15) 

^ {5'yl)U,s^, + {5<jl)\^^^„. (A16) 
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Let us notice that there is no mixed contribution since the cross-correlation {6(j)i£,a) = 0. Using the correlation 
function of the modulus noise, the term (^fi) can be further simplified, namely [76] 



Let us now turn to the second order correction. The solution for S<T2 can be written as [t^ 



6ff^((/>ci) 3i7-^((/)ci)J 3V;;(crci)ff^((pci) 



-(5o-i(t) + — ,,,, . 50i(t)Co-Mt. (A18) 



As expected, one sees that ba-i is quadratic in the noises. 

From the above expression, one deduces that the mean value of ba2 is non-vanishing and is the sum of various 
terms [ill 

where the last term in Eq. (jAlSP does not contribute because {S(l>i£,a) = 0. Had we not taken into account the 
stochastic behavior of the inflaton, only the term {Sa2) \scrf(^„) would have contributed. 

At this stage it is interesting to compare the previous considerations to Ref . ji] . In particular the above approach 
is a perturbative one and a relevant question is its domain of validity [t^, [ll]. In fact, as shown in Ref. [§], it turns 
out that, for the potential m^cr^/2, the Langevin equation can be integrated exactly. The solution reads 

a{t) ^ aine^ ^4. "^/(3H)dr ^ 1 ^ ^3/2 (^^^^^^^^j; ™^/(3//)d,^^ (A20) 
From this expression, it is easy to compute the variance. One obtains 

This is exactly the result obtained in Eq. (jA17|) . Therefore, although perturbative in nature, the approach used before 
has in fact a wider domain of validity and, in the specific case treated above, can be used even if the corrections are 
not small. The perturbative approach is also interesting for two reasons: firstly, it allows us to take into account the 
stochastic behavior of the inflaton field (even if, most of the time, we will show that these corrections are negligible). 
Secondly, the method of Eq. (jA17p rests on one's ability to solve exactly the Langevin equation which is possible only 
for a quadratic potential for a. If the potential is different (that is to say, not quadratic), only the method used here 
allows us to derive explicit results. 

We now turn to the calculation of the various corrections presented above in the following specific inflaton and 
modulus potentials. 



1. Chaotic rn^cj)^ inflation 

If we assume that = m^0^/2, then, in the slow-roll approximation, the classical evolution of the inflaton field is 
given by the following expression 




where N is the number of e-folds defined by iV = ln(a/ain), ajn being the initial value of the scale factor at the 
beginning of inflation, and 4'ci{N = 0) = 0in. The model remains under control only if the energy density is below the 
Planck energy density. This amounts to the following constraint on the initial conditions 0i„/Mpi < 87r\/2Mpi/m0. 
Inflation stops when the slow-roll parameter ei — —H/H^ is equal to unity corresponding to (/)cnd — ^/2M-p\. As a 
consequence, one can easily check that the argument of the square root in Eq. (|A22p remains always positive. The 
total number of e-folds during inflation is simply given by = / MpxY / A — 1/2. This number can be huge if the 
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initial energy density of the infiaton field is close to the Planck energy density. Finally, the inflaton mass is fixed by 
the WMAP normalization 

Qrms-PS ^ 1 ^ 1 ^4'* 

where the cosmic microwave background quadrupole is given by Qms-vs/T ~ 6 x 10^^ and a star denotes the time 
at which the scales of astrophysical interest today crossed the Hubble radius during infiation. Using the fact that 
ei* = 1/{2N^, + 1), where TV* is the number of e-folds between the time at which the physical scales left the Hubble 
radius during inflation and the end of infiation, one obtains 

127^^/T0 Qrms-PS r>-6 



Mpi 2iV* + 1 T 



7 X 10"" , (A24) 



where we have used N^, ~ 50. This also implies the Hubble constant at the end of inflation, i?i„f — m^/\^ ~ 
10" GeV. 

Using the perturbative presented before, one can solve Eq. (jAl[) (let us notice that the Langevin equation can 
be solved exactly only in the case of a quartic potential) and determine the quantum behavior of the inflaton field. 
Through straightforward albeit lengthy calculations, one obtains (ttI. Pzsj 

Note that the classical trajectory obeys (p < during inflation, so that the above is positive as it should be. 
Similarly, the correction to the mean value reads [73, \78^ 

9 

W.)--Y^^(0eV^?„). (A26) 

As before, since 4> < 0in during infiation, the quantity ((5(/>2) is positive. 

Concerning the modulus, we will consider two possible potentials, in the spirit of previous sections: one with fixed 
mass Va = m2(T2/2, and one with typical supergravity corrections. 



a. Modulus potential: Vo- ~ 



Let us first determine the classical trajectory of the modulus during infiation. Solving Eq. (IA2p without the noise 
term leads to the following solution 



(A27) 



In practice, one has 3> and, therefore, the modulus evolves slowly during inflation. Let us also notice that 



mj^/m-^ ~ 2/3{m„/Hinf^in) Nj, and, therefore, the limit rria/r 



corresponds to ma- ^ Hini, m-^T Even 



though the modulus is completely frozen in the limit nia/'m^ —> 0, it is necessary to take into account its evolution 
when computing the integrals in Eqs. (jA16|l . (jAl8|l . since its displacement can be non negligible if infiation lasts long 
enough. One then finds: 



ml 



192^2 A^4^ ^ _ 



^1/ (2m|) 

which, in the limit rria /m^ — > and (pd <^ (pin reduces to: 



(Sal) 



1927r2M4j 



4 4-2m^/mJ 2m^/m^ 



8^2 T 



(A28) 



(A29) 



Note that Hin refers to the Hubble constant at the onset of inflation, which differs (by a/2A^^) from the Hubble 
constant at the end of inflation, noted -ffi„f and used in the rest of our analysis. The two formulas (IA28P and (|A29p 
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are identical to Eqs. (23) and (24) of Ref. 0. The interpretation of Eq. (|A29[) is as follows. If the field is light, then 
the classical drift in the Langevin equation can be ignored. Then, it is easy to show that the Langevin equation can 
be integrated exactly and, as a consequence, that (cr^) = / H^dN/ (47r^). If the Hubble parameter is approximatively 
constant, then the previous expression reduces to <^cr^) — / (Att^) . However, as well known in the case of chaotic 

inflation, the Hubble parameter can change significantly during inflation. Therefore, the integral has to be evaluated 
exactly. When this is done, this produces the additional factor 1/2 present in Eq. (jA29p . 
Let us now turn to the second contribution originating from the inflaton noise. It reads: 



/A 2\ m^ggi (0ci-0in)'(</'ci + 3<^in) ^, ^Ti^^fn 1 mj a^, 



where the last expression is valid at the end of inflation. As a consequence, we note that {^'^i) ^^^g^^ is negligible in 
comparison with (Saf)^ , since ^ and ad/Mpi <C 1. 

Let us now calculate the correction to the mean value of the modulus field. The various contributions Sa2 amount 
to: 



miaci / 1 ,2 3 



192^2^/4^ \^ 2*^"' 2^'" 0ci 

™^f^cl /I 2 , 1 



{6a2)\sal(i.) = 0, (A31) 

the last result being obtained because V^" = in our case. Using again the fact that, as inflation proceeds, 0ci ^ 0in, 
one finally obtains 

Noticing that {5(72} is maximal at the end of inflation, i.e. for (f)c\ — \/2Mp\^ one can easily demonstrate that the first 
order correction is always the dominant one. 

Therefore, using the value of the inflaton mass obtained from the WMAP normalization, the expression of the initial 
value of the inflaton field in terms of the total number of e-folds, and setting > 60, one obtains the following 
lower bound: 

(Scrfy^'^ > 3.9 X 10-5 Afpi. (A33) 

The above considerations are valid provided one does not enter the regime of eternal infiation where the pertur- 
bative treatment of the quantum behavior of the infiaton field breaks down. Eternal inflation starts if (jjin ^ 
(24)^/''(to^/Mpi)-^/^. Therefore, the above calculations are applicable if < 2.8 x 10^ which implies that 
^/{Saf)/Mpi < 1.8. As was studied in this paper, such values of a are anyway excluded since the amount of en- 
tropy perturbations is too large to be compatible with the cosmic microwave background data. 

Finally, it is also interesting to investigate what happens if one considers a more complicated potential for the 
modulus. As mentioned before, the perturbative approach used here allows us to determine the stochastic behavior 
of a even if the potential is not quadratic. Therefore, let us consider the case where 

V.W^^^. (A34) 

where A„ is a dimensionless constant. The integration of the classical equation of motion leads to the following solution 



^ " K Mpj / (Pel 



+ (2-n) , " ^' In ^ 



1/(2-") 

(A35) 
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In the limit A„ — s- 0, the modulus field is almost frozen. Then, one can now compute the variance due to the modulus 
noise. According to Eq. (|A17p . it reads 



87r2 ^ V^Pi 
4- 3n 



2-2n 



4(3"-4)/(2-n) 



4(n - 1)!"^^ 
A„(2-n)M, 



■ exp 



PI 



4(n - l)!rn| 



2-n 



X S7 



2-n 



PI 



2-n 



4-3n 



2-n 



,-4 



fcl 



PI 



2-n 



(A36) 



where 7(a,x) = /q" *di is the incomplete gamma function. 



6. Modulus potential: Va = c^H^a^ /2 

In this subsection, one considers the case where the modulus mass receives supergravity corrections of the form 
cfH^i assuming q < 3/4, during inflation. A comment is in order at this point. One uses the notation a in order 
to emphasize the fact that the supergravity corrections to the modulus potential are not necessarily the same during 
inflation and during the post-inflationary epoch. As a consequence, one expects Ci ^ c. This is of course the same for 
the minimum of the potential which, during inflation, is not necessarily equal to (Jq, the minimum of the potential in 
the post-inflationary epoch. In the rest of this appendix, we work in terms of the held displacement (with respect to 
the inflationary minimum) rather than in terms of the fleld itself. 

Then, if one assumes that the supergravity corrections dominate, then the modulus classical motion reads: 



Cci = cTin exp 



exp 



3Af2j 



dip 



(A37) 



and for micf)^ inflation, this gives: 



(Tin exp 



12M2j 



(A38) 



Very quickly the argument of the exponential becomes — cfcpf^/ (12 Mp^) ~ — c? A'^/3 



(m^/i/i„f,i„)2iV^/3 and. 



therefore, as previously, one recovers that the massless condition is given by rricr ^ Hinf^nN^^^"^ . Since > 60, this 
condition is now always violated (at least provided that Cj > 0.22). 

As mentioned above, Eqs. (IA14p . (jAlSp need to be corrected to account for the explicit time dependence of the 
potential, which enters through the H prefactor. Making this dependence more explicit in Eq. (|A13I) . one solves these 
equations as 



Sai = -fa.i dTH'iT)6Mr) 



dr exp 

27r 



drj H(t]) 



(A39) 



and hence 



167r2M2jcf 



6M|i 

cl 



cf(j)l 



(A40) 



Therefore, one obtains the Bunch-Davies result (with the mass cfH^) corrected by a factor the value of which at the 
end of inflation is l-|-3/c?. This result agrees with that obtained in Ref. 9]. In the limit q 0, one of course recovers 
the result Eq. (|J29|) . 

The contribution to the variance due to the inflaton noise is given by 



'-i'^cl"V 

' 34567r2M|i 



1 

30' 



5 



(A41) 



This term is negligible compared to the previous one since the vev of the modulus is small (in Planck units). 
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Let us now turn to the corrections 802 to the mean value. As explained before, one needs to modify the results 
above to take into account the fact that the mass is explicitly time-dependent. Straightforward calculations lead to 



(^0-2)5^2 
This can be integrated to give: 

('^Cr2)5rf,2 = 



H"{t) H'^t) 



3 

. 



dr 



tin 



H{r) 



iJ'(r)drei- <^^h(^) ^^^^(^^^^^ ^^^^ 



11527r2M|i 
11527r2M|i 



69127r2M|i V 30 



= 



(A42) 



(A43) 



Given that fici/Mpi <C 1, it is easy to see that these contributions are subdominant. Therefore, the main contribution 
is the one given by Eq. (jA40|) . Expressed at the end of inflation (0ond = V^Mpi) and normalized to the cosmic 
microwave background, this expression leads to the following constraint on the value of the modulus: 



(Sal) 



1/2 



> 0.79 X 10" 



PI ■ 



(A44) 



2. Small field and hybrid inflation 



In this section, we turn to another type of inflationary model. We now consider a potential of the form: 

PI 



4> 



M 



1 



e I - 



(A45) 



with e = ±1. Such a potential gives rise to small field inflation if e = — 1 or hybrid inflation if e = +1. Since all 
integrals cannot be carried out exactly in this case, we provide the results to leading order in Out of simplicity, 
we use the notation $ = (/)ci/m ^-i^d similarly for <i>in. 

Let us now discuss small field inflation in more details. For this model, the slow-roll trajectory is only known 
implicitly. It can be expressed as 



N 



2p 



$2 



$2 



p-2 



2-p 



.$2-P 



P 



(A46) 



If p = 2 the singular terms must be replaced by a logarithm. From the above formula, one deduces that, given that 
(/"in ^ 0cnd, the total number of e- folds can be written as 



1 



Af li p{p - 2: 



2-p 



(A47) 



In this class of models, the end of inflation occurs by violation of the slow-roll conditions. If ii/Mp\ <C 1, it happens 
at 'i^end — [2/j,^/(p^Mp[)]^/*^^P-^^. In this regime, the two first slow-roll parameters are given by 



El* 



exp 



-47V, 



M, 



PI 



£2* 



(A48) 
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for p = 2, while, for p ^ 2, one has 

-2(p-l)/(p-2) 



El* 



p2 fMpi\ 2 



PI 



e2* = ^^. (A49) 

Jy^, p — 2 



If /i/ilfpi > 1, then the above formulas are no longer valid and the slow- roll parameters must be evaluated numerically. 
In particular, in the limit fi/Mpi 3> 1, one has 0cnd — A*- 

Our next step is to deduce an expression of the mass scale M. One obtains 

\ 7on_2„2'^rms-PS r ^„ oMvr l-2(p-l)/(p-2) / A* 



Mpi 



720.V^^ [Pb- 2)7V.]-^(^-i)/(^-^) (^^j , (A50) 



where we have assumed p ^2. This expression is valid only if jji/M^i ^ 1. If /i/Mpi > 1, one can show, see Ref. [79| . 
that the end of inflation occurs when the vacuum expectation value of the inflaton field is of the order of fi and, as 
a consequence, that (M/Mpi)* ~ O (Qrms-ps/^^) ■ ^^^^ notice that, contrary to the case of chaotic inflation, 

the Hubble parameter is approximatively constant during inflation and given by iJinf ~ MV(y3A/pi). Therefore, if 
M ^ ^Pi, -ffinf — 10^^ GeV, but TJinf can be much lower if /Lt < Mpi. 

Let us now discuss the constraints on the free parameters p and /i. As shown in Ref. [79], there is no prior 
independent constraint on the index p. However, if one adopts the theoretical prejudice that the vacuum expectation 
value of the inflaton field must be smaller than the Planck mass, then must be smaller than Mpi and then one 
can demonstrate that the case p = 2 is slightly disfavored by the cosmic microwave background data. In small 
field inflation, the energy scale of inflation can be very low. But, it can not be smaller than, say, the TeV scale, 
T4/mp[ > lO^®"'. If p = 2, this implies that ^/Mp\ > 1.25. In this case, the formulas (|A48|) are no longer valid but 
one can show that the spectral index is still compatible with the data in the regime /i/Mpi > 1 or even /z/Afpi ^ 1. 
If p =/: 2, the parameter fi is globally unconstrained. The constraint mentioned before implies 



M . j4x 10-64 /Q,,„3_ps 



_2 ^ (p-2)/(2p) 



>}lJL^n±I^\ [p(p-2)7V,f(^-^)/(P-H . (A51) 



Mpi ^ 1 45p2 \^ T 

This leads to fi/Mpi > 0.16 for p = 2.1, ^/A/pi > 3.4 x 10"^ for p = 2.5, /^/Mpi > 1.5 x 10"*^ for p = 3, /i/A/pi > 
9.9 X 10^^^ for p = 4 etc .... In this case, one sees that /z/Mpi can be small. For p 2, Eqs. (|A49p indicate 
that £2* no longer depends on fi. The spectral index remains compatible with the data in this regime because 
£2* ^ 0.04(p— l)/{p— 2) still lies in the 2(t contour whatever the value of p. As already mentioned, if, on the contrary, 
/i/Mpi ^ 1, Eqs. (|A49p are no longer valid but one can also show that the model is still in agreement with the data. 
The conclusion is that, provided the energy scale of inflation is above the TeV scale, the two regimes fi/AIpi ^ 1 and 
fi/Mpi > 1 are still compatible with the data and, hence, the parameter /i remains basically unconstrained [79j . 

Let us now briefly discuss hybrid inflation. A crucial difference with small inflation is that inflation no longer stops 
by violation of the slow-roll conditions but by instability. This means that there is one more additional parameter 
namely the value of the inflaton field at which the instability occurs. This makes the analysis more complicated since 
the parameter space is enlarged. For this reason, although we give all the necessary expressions, we have chosen in 
this paper to skip a detailed investigation of this case. 

Straightforward calculations lead to the following results for the variance and the mean value of the inflaton field 

(^^?) - - ZtC (*^-^ - a>f-^^$^-^) . (A52) 



Similarly, 



127r2Af4j ep(4 - 3p) 



p-2.^ A~3p 



247r2A/4j ep{A - Sp) 



+ Z 2il$P-l$2-2p ^ _ l)$2p-3^4 



4-3p 



(A53) 



These results agree with those of Ref. [77| to leading order. 

As regards to the modulus, we again consider two possible potentials. We consider these two possibilities in the 
following. 
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a. Modulus potential: Va ~ rn^a'^ /2 
The classical trajectory of the modulus reads: 



o-ci = o-inexp 



9 9 



ep{2 - p)M^ 



(A54) 



Let us notice that the argument of the exponential is always negative. We now discuss the case of small field inflation. 
The argument of the exponential in Eq. (|A54p is always dominated by the following term 



p(p-2)M4 - 3 V^finf 

The argument of the exponential can be small or large depending on the total number of e-folds during inflation and 
on whether the modulus field is light or heavy. Using the WMAP normalization, one can also express the argument 
of the exponential in terms of the parameters of the model. In the regime fi/Mpi <C 1, one obtains 

Let us give a few examples for different values of p with ~ lO^GeV and N^, — 50. For p ~ 2.1, this gives 
i3.78N^{^i/Mpi)-^^, for p = 2.5, one has 6.29 x IQ-^^N^ifi/Mpi)-'^" , for p = 3, one obtains 3.71 x IQ-^^N^ifi/Mpi)-^ 
and for p — 4, this leads to 2.64 x 10~^^A'^ (/i/A/pi)^'*. Let us be more precise and give some numbers. If p = 2.5 
one can easily check that the argument of the exponential is always large. But this conclusion can be modified if one 
considers other values of the parameters. For instance, if p = 3 and 1.5 x 10^^ < ^/Mp\ < 8.4 x lO^^N^^^, then 
the argument of the exponential is large but, if 8.4 x 10~^N^^^ < /i/Mpi <C 1, it is no longer the case. Let us also 
notice that the previous considerations are valid only for values of such that the number 8.4 x 10~^A'^/^ remains 
small otherwise the formulas used here would not be valid. This means A^^ <C 1.84 x 10^ which is not so restrictive. 
The same analysis is true for p = 4, the corresponding intervals being 9.9 x 10^^^ < ^/Mp\ < 4 x 10~'*A''^/^ and 
4 X IQ^'^N^/'^ < ^i/Mpi <^ 1. In the present case, the total number of e-folds must satisfy <€. 3.9 x 10^. On the 
other hand, in the regime where /x/Afpi > 1, the argument of the exponential can now be written as 

and is small, unless we choose very large values of A^^ . 

We are now in a position where one can compute the variance of the stochastic motion using Eq. (jA16p . This leads 
to the following result 



ep{2 - p)M^ 



(A58) 



Let us discuss this result. If the argument of the exponential in Eq. (|A54[) is large (in absolute value), then the 
exponential becomes negligible and 



("^^i)?.?. - 247r2M4,m2 " 8^2;;^2 ' (^^9) 



where we used Hf^^ ~ M^/(3Mpj). It was suggested in Ref. [9] that this formula describes the small field case. But, 
as was noticed before, the argument of the exponential in Eq. (|A54p can also be small and, in this case, one has 



We notice that this last expression is similar to Eq. (jA29p . The only difference is that we do not have the presence of an 
additional factor 1/2 with respect to the standard case. This is because, in the case of small field inflation, the Hubble 
parameter is almost constant and, hence, one obtains the de Sitter result. We conclude that if rria ^ H^f / \/ , then 
the variance is given by Eq. (|A60p while if Hinf/\/N^ ^ ra^ < -ffinf, it is given by the Bunch-Davis expression (|A59p . 
Finally, if m„ > -ffinf , then the variance becomes negligible. 
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As was done for the case of chaotic inflation, one can also express the variance in terms of the parameters of the 
model. Let us start with the case fi <^ Afpj. In the regime where the Bmich-Davis term dominates, the variance can 
be written as 



iUpi 



1.63xl0V|p(.-W*-"'"-'(^)" (jI;) . (AM) 



For p = 2.5 and m„ = lO^GeV, this gives (Saj)^ ^ /M^^ ~ 1.79 x 10-^{^/Mpi)^^ which implies that 7.6 x 10"^^ < 
(Scrl)^ ^ /Af|j < 1.79 X 10^29^ On the other hand, if p = 3 and m„ ^ lO^GeV, this leads to (Saf)^ ^ /M^j ~ 5.15 x 
IQ-^ill/MpiY^ which means that 6.7x 10^^^ < {Saf}^ ^ /M^j < 6.35x IQ-^Oivs. Finally, if p = 4 and = lO'^GeV, 
this means that (Saf)^ ^ /M^^ ~ 0.01(^/Mpi)8 which implies that 9.2 x 10"^^ < {5(^1)^ ^ /M^^ < 6.5 x IQ-^OiV^. 
We see that we always find a very small contribution. It is easy to see that it originates from the fact that the energy 
scale of inflation is very small and that we always have M^/Mpj <C mcr/Mp\. Let us now consider the other regime, 
given by Eq. (jA60p and where the Bunch-Davis term is subdominant. In this case, the variance can be expressed as 

(Sa^) / \ 2p/(p-2) 

^ 2.16X 10-X/[p(p-2)iV*]-''^-''/(^-') -r^) . (A62) 

Therefore, if to^ = lO^GeV and p = 3, one obtains (Scrl)^ ^ /M^^ ~ 3.84 x lO^^^TV^ (^i/Mpi)*^. This gives the 

following range of values: 1.35 x lO^^^N^^ < (Sal)^ ^ l^'lpy S5 3.8 x IO^^s^y^ This leads to a significant constraint 

only if the number of e-folds is large (remembering that, for p — 2>, the formulas applies only if ^ 1.84 x 10^, see 

above). Taking this upper limit, one arrives at ^JJSaf}^~^ / Mpi ~ 8.3 x 10~^ which, besides being an extreme case, 

remains small in comparison to what was found in the case of large field models. If p = 4, one has (^Saf)^ ^ l^p\ — 

5.4 X N^{ii / Mp\Y and the same conclusion can be reached. However, the extreme case evoked before becomes 

more significant as one obtains \J {5aD^ ^ /Mp\ ~ 1.4 x 10^'"'. 

Let us also consider the regime where ^/Mp\ > 1 or /i/Mpi 3> 1. As mentioned above, in this situation, the 
argument of the exponential in Eq. (|A54p is always small and iJjnf ^ 10^^ GeV. This means that one is always in the 
case where the Bunch-Davies term is subdominant. Then, given that > 60, one obtains the robust lower limit: 

{^<^D]'l > 4.2 X 10-*^ Mpi . (A63) 

Let us now turn to the other contribution to the variance originating from the fact that the inflaton noise. It can 
be expressed as 
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The analysis is complicated by the fact that the dominant term depends on whether p < 4 or p > 4. Since, previously, 
we have mainly considered situations where p < 4, we will restrict ourselves to this case. In addition, very large values 
of p appears rather unnatural form a high energy physics point of view. Therefore, in the case of small field inflation, 
the dominant term can be written as 

/A^2\ _ /i-^m>gi 1 

- 247r2M4iMX4-3p)(4-2p)(4-p) ' ^ ' 

In the regime where ^/Mpi <^ 1, the term can can be re-expressed as 
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If p = 2.5 and = lO^GeV, then one obtains (SctI) g^^g^^ /M|j ~ 1.9 x 10-*^{aci/Mpi)'^{ii/Mpi)-^ . But we have 
seen before that, Hp = 2.5, the argument of the exponential in Eq. (|A54[) is always large (and negative). This means 
that the modulus is exponentially killed during inflation and, hence, negligible at the end of inflation. In this case, the 
contribution {5<jf)g^ A^pi is also negligible. If p = 3, the story is slightly different. Indeed, for rria — lO^GeV, one 

obtains (^cr^)^^^^^^ - 6 x 10-^\aci/Mpi)^{n/Mpi)-^/^. If 1.5 x 10"^ < /i/A/pi < 8.4 x lO^^N^/'^, the argument 

of the exponential in Eq. (|A54[) is large and the same conclusion as before applies. But when 8.4 x lO^^N^/^ < 
fi/Mp\ <C 1, the argument is small and the modulus is almost frozen, CTci — (Jin- In this case, this implies that 
1.9 X lO-^^{ain/Mpif < {Scrf)g^^g^^ < 7.8 x lO-'^^Nl/^{anjMpi)^ . We see that this contribution remains very 

small. 

Let us now consider the case where fi/Mpi > 1. As already mentioned, this means that (/)cnd — fJ-- In addition, we 
have shown before that, in this regime, the modulus is almost frozen during inflation. As a consequence, the variance 
due to the inflaton noise can be expressed as 
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We conclude that this contribution is negligible unless one takes very large and unrealistic values of /x. 

After having estimated the variance, one can now calculate the correction to the mean value. The non-zero 
contributions to {5a2) read: 
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We notice that the corrections to the mean value are of the same order of magnitude as {Scrl) ■ Therefore, 
following the above analysis, one can safely conclude that these contributions are negligible. 



h. Modulus potential: = 



If the modulus mass is dominated by a Hubble term contribution from supergravity effects, one obtains the following 
results for the classical trajectory and the quantum corrections: 
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At this stage, one can reproduce the discussion of the previous section and study the argument of the exponential. It 
is easy to show that it can be expressed as 
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We see that the above formula is nothing but Eq. (|A54p with the time dependent mass CiH. The argument of the 
exponential can be large or small depending on the parameter Ci and the total number of e-folds. But, as already 
mentioned in the section devoted to chaotic inflation, if Ci > 0.22, then it is always greater than one given the fact 
that iV^ > 60. 

The variance of the first order correction is made of the combination of the following two terms, as above. The 
term due to the modulus noise reads 
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If the argument of the exponential in Eq. (jA71|) is large, then one has 

and one recovers the Bunch-Davis term with a mass cHi^f. On the contrary, if the argument of the exponential in 
Eq. (|A71[) is small, then the expression of the variance reads 

and one recovers the de Sitter result. In particular, the term cf has cancelled out, as, in the corresponding situation, 
the modulus mass nia cancelled out in the previous section. 

As before, one can also express the above results directly in terms of the relevant parameters. Let us start with the 
situation where the argument of the exponential is large. In the regime /z/Mpi <g; 1, this gives 
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If p = 2.5, this gives ((5cr?)^ ^ /Af^j ~ 3.39 x IQ-^^ct'^ (^/Mpi)^° and if p = 3, one has {Saf) ^ ^ /M^j ~ 5.76 x 

lO^^^Cj"^ {fi/AIp\)^. Given that the above estimates are valid in the regime where p,/M-p\ is small, we conclude that 
the variance is always very small unless Ci takes tiny and unrealistic values. 
In the regime /i/A/pi > 1, one obtains: 

(5al)]'l 6.8 X 10-^^ . (A77) 

Contrary to the previous case, the above formula indicates that the quantum effects can now be significant, especially 
when the parameter q is small (q < 0.22). 

If the argument of the exponential in Eq. (|A7ip is small, then it is easy to see that the variance is now equal to 
2iV.j,cf/3 times the variance given in Eq. (jA74|) . In the regime fJ,/Mpi <C 1, the factor N^c? is unlikely to compensate 
the smallness of the variance obtained in Eq. (IA74P (unless the total number of e- folds is huge), and we conclude that 
the quantum effects can become arbitrarily small in this regime. 

On the other hand, if n/Mp\ > 1, then one obtains (Saf)^^^^ /M|i ~ 3 x IQ-^^N^ and, given the fact that the 
total number of e- folds must be larger than 60, one obtains the following lower bound 

(S<^f)yi > 4.2 X 10"^ i\/pi . (A78) 
Let us now consider the second contribution to the variance due to the inflaton noise. Its expression can be written 

as 
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As discussed in the previous section in a similar context, the amplitude of this term depends on whether p < 4 or 
p > A. Here we restrict ourselves to the case p < 4. Moreover, if the modulus significantly evolves during inflation, 
then the above contribution becomes negligible. So we assume that cfiV^/3 ^ 1 which implies that a is small. In 
this situation, if /z/A/pi <C 1, one has 
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This term is obviously tiny, in particular because the power of the term /i/Mpi is positive. Hence, the contribution 
due to the inflaton noise can be neglected. 

The case where /z/Afpi > 1 remains to be studied. Straightforward considerations lead to 
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This contribution is small unless /i/A/pi ^ 1. 

Finally, the (non-zero) second order corrections to the mean value of the modulus vacuum expectation value read: 
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It is clear that the order of magnitude of those terms is similar to {^<^i)g^ ■ Hence one concludes that the corrections 
to the mean value can be safely neglected, at least when fi/Mp\ is not too large. 



APPENDIX B: DETAILS OF THE SOLUTIONS TO THE EQUATIONS OF MOTION IN PRESENCE OF 

HUBBLE SCALE MASS CORRECTIONS 

In this appendix, we study in details the solutions of the modulus equation of motion in the post-inflationary epoch 
in the case where the potential receives supergravity corrections. 



1. Case +c^H'^/2 

The solution to the equation of motion has been obtained by Lyth & Stewart [l2|. The field equation can indeed 
be put in the following form: 



with the following definition: 
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A prime denotes derivative with respect to m^t^ while p characterizes the global equation of state of the Universe: 
H = p/t^ with p = 2/3 for matter domination or p = 1/2 for radiation domination. The above field equations takes 
the form of a Lommel equation. Out of convenience, we write the general solution as the sum of a homogeneous 
solution with constants ai and q;2, and a particular solution given in terms of a hypergeometric function: 
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The initial conditions are defined at t = p/Hinf as follows: a = dinf and cr' ~ 0. As a consequence, the coefficients 
Qfi and a2 can then be expressed as: 




ai = (fTinf - Co) 



"2 = (o-inf - O-q) 



Obviously, if dinf = ctq, meaning that the vev of a at the end of inflation corresponds with the high energy minimum, 
then ai — a2 — 0. The solution then reduces to the particular solution. 

We now study the behavior of the solutions according to the values of the parameters. Let us start with the situation 
where is real, namely c < (/i + At early times, mo-t <SC 1 (hence H ^ rua-) and at late times, nia-t ^ 1 

(hence H ^ mo-), the particular and homogeneous solutions evolve differently. At early times, the particular solution 
remains constant, 



(B7) 



Therefore, before the onset of oscillations, the particular solution does not undergo any redshift. On the contrary, at 
late times, one obtains 

apart(0 - t^o(m,t) — /^rf- + ^^jrf-- + ^l^jsinfTO,i+--^-^7rj , m.i » 1 , (B8) 



and, in this regime, one has that crpart(0 ^ a~^^^ regardless of p. Since, in the hmit m^t '3> 1, the supergravity 
corrections are by definition negligible, one can write — 'Ti^tTpjjj.j Qgj.i(a/aosci)~^/2. From the above formula, one 
deduces that the quantity (Tpart,osci can be expressed as 
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;/2 only. Indeed, 



Let us notice that the modulus energy density at ma-t ^ 1 (or H 

at nia-t ^ 1, the term c^H^ can not be neglected and participate to the energy density at the onset of oscillations. 
Therefore, strictly speaking, rn^a^^^^ osci/2 only represents the contribution of the term to^(T^/2 ai H — However 
what really matters is the modulus energy density at the time of dark matter freeze out. In this regime, the supergravity 
corrections are negligible. As a consequence, one only needs to take into account the contribution due to the term 
TO^(T^/2 and to express it in terms of its value at the onset of oscillations. Then the contribution of the particular 
solution to the energy density contained in the modulus is given by 
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Note that this result differs from that obtained for the purely quadratic potential: in the latter case, r^o-^osci is controlled 
by the value of a at the end of inflation, while here, it is controlled by the value of the local minimum ctq generated 
by supergravity corrections. 

Another issue concerns the values of the parameters p, and u that should be inserted in the above equation. If 
Tih > Tosci, then one should use /j,rd and i^rd, where quantities indexed with "rd" (resp. "md") are to be evaluated 
for a radiation dominated (resp. matter dominated) era, that is to say with p = 1/2 (resp. p = 2/3). The justiflcation 
is the following one. Just after inflation, the universe is matter dominated and the solution (|B3[) with p — should 
be used in order to describe the evolution of a. Since, in the present case, we study a situation where the onset of 
oscillations happens after reheating, one is in fact in the regime m„t <^ 1 during the whole matter dominated era and, 
as a result, the particular solution remains constant and equal to ctq during this phase. Then, after the reheating is 
completed, the universe becomes radiation dominated and, in principle, one should solve again the equation of motion 
but, this time, with p — 1/2 rather than p = 2/3 and with "new initial conditions" inherited from the reheating era. 
However, at the beginning of the radiation era, we are still before the onset of oscillations (i.e. still in the regime 
rria-t <^ 1) and the solution with p = 1/2 is still constant since this asymptotic behavior does not depend on p. 
Consequently, the matching of the two solutions is in fact trivial. This is means that, in the radiation dominated era, 
the relevant solution is nothing but the solution (jB3p with p = 1/2. Therefore, at late times, it is sufficient to use 
Eq. (jBSp with p = 1/2, hence the above claim. 
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If Tosci > Tlh, the situation is slightly more complicated. This times, the onset of oscillations occurs during the 
matter dominated phase and the matching between the solution with p — 2/3 and p — 1/2 should be done in the 
rricrt 3> 1 regime. However, as before, the two solutions are the same because a scales as a"'^/^ independently of p, 
see above. As a result, one can ignore the matching and still work with the solution with p = 2/3 at late times in the 
radiation dominated era. Therefore, if Tosci > Trh, one should use fiMD and vmd- 

Let us now consider the homogeneous solution. At early times, it can be approximated as: 
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where one has used the expressions (|B5p and (|B6p . In particular, the last terms in Eqs. (|B5p and (|B6p become 
negligible which implies ai ~ 02 and, therefore, the terms proportional to J^^ cancel out. One notices that the 
homogeneous solution does not remain constant prior to the onset of oscillations as it was the case for the particular 
solution. One finds that po-.kin and Pa.pot scale as cx i2(-^-2+i>) p^j, ^ matter dominated era, one therefore finds 

On the other hand, the late time evolution 
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As expected, one finds the standard evolution at late times, with a ^ since /i + 3/2 = 3p/2. In one uses the 

expressions of ai and a2 given in Eqs. (|B5p and (|B6p . this can be re- written as 
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As before, one should distinguish two different cases. If Tosci > Th, then the onset of oscillations occurs in the matter 
dominated era. As a consequence, at the end of the reheating stage (in the regime mat 3> 1), the homogeneous 
solution is given by 
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Then, one should match this solution to the solution valid in the radiation dominated era. However, as in the case 
of the particular solution, since one is in the regime mat 3> 1, the two solutions scale as a"'^/^ and the matching 
becomes trivial. As a consequence, one finds 
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Finally, the case Tih > Tosci remains to be treated. In this situation, one needs to take into account the change in the 
values of /i and v and there is no other choice than performing the matching explicitly at the time of reheating. This 
time, the matching turns out to be non trivial. The whole matter dominated era occurs in the regime mat <^ 1 and, 
to leading order [using again the expressions (jBSP and (|B6p ] one obtains that 
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In the radiation dominated era, the solution is given by = (3i{mat)^^^'"°^^ Ji,^^{mat) + (32{mat)^^^^°^^ J^i,p^^{mat) 
and one should determine the coefficients /3i and /?2 by matching, at reheating, the above solution to the values (jB16p 
and (|B17p . Straightforward calculations show that pi ^ 132- Then, using the asymptotic behavior of the Bessel 
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function J^no i'^a't) that appears in the dominant branch of the solution in the regime m^t ^ 1, i.e. after the onset of 
oscillations, one can check that this solution scales as a~^/^, the proportionality coefficient being, as before, directly 
related to ri(T,hom,osci- One finds 
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As expected, one can check that, in the limit where all the quantities labelled "RD" becomes equal to their counterparts 
labelled "MD" , the above expression exactly reduces to Eq. (jBlSP . 
Summarizing our result, one has 
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Let us emphasize that, in order to obtain the above expressions, we have assumed that the total energy density at the 
onset of oscillations is only made of two pieces, one originating from the particular solution and the other from the 
homogeneous solution. This means that we have neglected the cross terms that, in principle, should contribute. This 
is justified by the fact that, in practice, one solution always dominates the other. Which one is dominating depends 
on the region explored in the parameter space. It is also convenient to define 

A, . ^^^r (^"" + ;"" + ^) r (^^"-;^° + ^) , (B21) 

and A2 expressed in the same way as Ai above, but with /imd and /iMD replaced by /ird and z^'rd; The coefficients 
Bi and B2 are given by 
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'^RD V ^RD / 

With these definitions, the above expressions (jB19p and (jB20[) exactly match the ones used in the main text, see 
Eqs. dMl) and (HOD. 

The other possibility is when v is imaginary namely c > (/i + For p — 2/3, this case corresponds to c > 3/4. 

Writing v = zi>, one has v > ^ growing with c. This brings in non-trivial modifications for the scaling of the modulus 
energy density. Following the same reasoning as before, one finds that the contribution due to the particular solution 
can be expressed as 
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As before, in the above equation, one should use quantities labelled "MD" or "RD" according to Tosci > T^h or 
^osci < Trh- The difference with Eq. (jBlOp comes from the fact that the numerical prefactors, which are of order unity 
if V is real, may become quite small if v is pure imaginary and its modulus is large. This point is elaborated further 
in more detail in the main text. 
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Turning to the homogeneous solution, we finds that it scales as follows. At early times, one no longer finds ai ~ a2 
and this his implies that 



a{t) ~ t 
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Keeping in mind that 1/ is imaginary, this gives the following scaling of the kinetic energy density: Pkin,hom oc ^^(-A'-a)^ 
and it is easy to see that at early times, the potential energy density scales similarly. At late times, the homogeneous 
solution evolves as 
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As expected and as it was the case before for the case of v real, one can check that a scales as oT^I"^ regardless of 
p. At this point, as was done before in the case where v was real, one has to distinguish the case where the onset of 
oscillations occurs before or after the reheating. Let us start with the case Tosci > ?rh. Since the oscillations start in 
the matter dominated era, one is interested in the regime m^-i S> 1 with p — Iji. The corresponding solution is given 
by Eq. (jB26p . Expressing the coefhcients a\ and 0L2 using Eqs. (jBSP and (jB6|l . one obtains the following result 
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where Tmd = ^'md ln[pma-/(2i?inf )] — 6md — ^'md, Bmd being the phase of the complex number r(l + ij>MD) and 
^MD the phase of 1 — i(/xmd + 1)/^'md- As usual, this function scales as a~^^^ . Then one should match this solution 
with the solution in the radiation dominated era (still in the regime mat ^ 1). However, the solution in the radiation 
dominated era also scales as a"'^/^. As a consequence, the matching is trivial. Therefore, using the above equation in 
the limit where ji'l 3> 1 (which simply amounts to replace, in the above expression, sinhx and coshcc by e^/2), one 
obtains 
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where we have ignored the remaining trigonometric functions since they give contributions of order one. It is remark- 
able that all the exponential dependence in has cancelled out in this expression. 

When Trh > Tosci, the onset of oscillations occurs in the radiation dominated era and one additional matching is 
required as explained previously. At reheating, we have 



t 



CT|rh 

do; 

rh 



(finf — CTq) 
— (Cinf — O-q) 



inf 



Hini 



1+ ( ^MD + 1 



cos 



inf 



-MMD— 2 



sm 



Hrh 

i>MD In 



+ «'md 

ffinf 



(B29) 
(B30) 



These equations should be compare with Eqs. (|B16|1 and (|B17|1 . As before, in the radiation dominated era, the solution 
is given by tr = Pi{mat)~^^°~^ J„ji^{mat) + I32(mat)^'^'^°^^ J-^-p^-^lmat), where t'RD is now a complex number. The 
coefficients Pi and /32 are determined by matching the previous solution to the values (|B29p and (jB30p . Then, the 
solution is completely specified. In the regime mat 3> 1, one obtains 
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As expected, this solution scales as a~^/^. Again, it is interesting to observe how any exponential dependence in v^X) 
is in fact exactly cancels out. It is also interesting to study the case where all the quantities labelled "RD" becomes 
equal to their counterparts labelled "MD" . In this case, lengthy calculations using trigonometric identities show that 
the term in the curled bracket reduces to 



cos *MD 



COS Tmd cos (^m^t ~ + ^md sin ^m^i — — ^ 



(B32) 



and recalling the definition of 'I'md (which implies that cos exactly cancels the square root in the above formula), 
it is easy to demonstrate that Eq. (jB3ip exactly reduces to Eq. (|B27p . Finally, one obtains 
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This completes our analysis of the case where v is complex. 

As before, the total llo-,osci is roughly given by the sum of rio-^hom.osci and rio-,part,osci- Summarizing our results one 
obtains 
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-} , r,h > ro.diB35) 



These equations are the counterparts of Eqs. and ([50)1 in the case where the quantity i' is complex. As was done 
before, it is convenient to define a prefactor ^3 (the equivalent of Ai, see before) by 



3 . i'MD 
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As before, if Tosci < Trh, one should rather use a prefactor Aa defined as ^3 above, but with all ^md and t'MD replaced 
by their values for the radiation dominated era, /xrd and ^rd- Finally, for the homogeneous solution, one introduces 
the coefficients 
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With definition The above expressions (|B34|1 and (jB35P reduces exactly to the expressions used in the main text, 
namely Eqs. and 



2. Case 

As explained in the main text, the potential can be written under the approximate form: 



A„ (a - fT„)4+« 



(n + 4)! 



(B39) 



PI 



with: 
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The vev tT„ indicates the (time-dependent) local minimum of the effective potential. 
As discussed in the main text, the potential can be further approximated by V{a) 



A„a"+V[(n + 4)!M, 



piJ 



the 



limit a ^ an- Then one can show, as follows, that the field will evolve to values of order (T„ in a fraction of an e-fold 
of order (tT/tT„)~^"+^)/^. In the following, we write (T„^inf the value a„ {H — i?inf)- Consider the field to be initially at 
rest, (Tinf = 0. Then, neglecting the damping term H& in the equation of motion, one can rewrite this latter equation 

as 
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CTinf 
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which shows that the field initially evolves on the e-folding scale announced above. As soon as the field starts to 
evolve, it converts its potential energy into kinetic energy; since this occurs in a fraction of an e-fold, expansion 
damping can be ignored and the conversion is nearly fully efficient: Pa,kin — ^(cinf) — ^(o')- Once the kinetic energy 
dominates over the potential energy, Pa-.Mn oc a~^, and one finds the approximate solution: 
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3-ffinf \ Hinf , 

assuming matter domination after inflation. The factor ^y2V{ainf) is introduced to approximate the value a at the 
time at which the kinetic energy of the modulus dominates over the potential energy. The ratio ^yV(a\^ / Hinf ~ 
Cinf (o'inf/o'n^inf)'-"^^-'^^, therefore the above solution shows that a evolves down to (T„ on the aforementioned e-folding 
scale. 

When (T ^ cr„ (and _ff ^ m„lc]. one can approximate the potential (jB39[) with the low order term +c?H^{a—a„\^ . 
In this case, the solution of the equation of motion for a is the sum of a particular and an homogeneous solution, 
namely 



ait) = a„cr„(i) -I- CThom(0 , 



with the definition 
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The homogeneous solutions can be expressed as follows. Reintroducing the notation: — (/i + 1)^ — p^(? , the 
solutions read differently according to whether v is real or imaginary. Let us start with the case v real. In this case, 
the solution normalized to negligible initial kinetic energy and initial value of the modulus ui reads (4^ : 
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The time dependence is contained in the dependence of H{t). Of course ai must be chosen so that the sum of 
Chom + fpart takcs the right value at the end of inflation. The question of the initial conditions in this case is discussed 
in more detail in the main text, see Sec. IIII B 21 Depending on the value of n, either the particular solution or the 
homogeneous solutions may dominate at late times. In general, however, one should expect the homogeneous solution 
to decay less rapidly, since this corresponds to: n < 2(— yu -\- v) / [n + l — v), and during matter domination, p = —1/2. 
Let us now study the case where v is imaginary. As usual, we write v = iv and, then, the solution reads p3] : 
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The homogeneous solution scales less rapidly with time than the particular solution if n < — 2/x/ (/x + 1). 



APPENDIX C: SOME CONCRETE MODELS OF MODULI EVOLUTION 



In this appendix, we present specific supergravity examples of inflationary scenarios which give different predic- 
tions of modulus evolution. Since the moduli fields are scalar fields with a flat potential which eventually get a 
supersymmetry-breaking mass of order of the gravitino mass the simplest example is the supersymmetry- 

breaking field or the Polonyi field itself (denoted S in what follows). Therefore, we consider a system consisting of an 
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inflaton sector and a separate Polonyi sector. In the rest of this appendix, we assume that the Polonyi sector is charac- 
terized by a minimal Kahler potential, Kg = jS'p, and a superpotential given by Wg = M| {f3 + S). Here /3 is a con- 
stant chosen as /3 = (2 — \/3) Mp\ so that the potential energy vanishes at its minimum S'min = {V^ — l) Afpi when the 

inflaton sector is absent. At the potential minimum the gravitino mass reads m^/2Mp\ — (^e^^min^'^Ws^i^) = M|e^^^. 
Below we mostly take Mpi ~ 1 but recover it when appropriate. 



Small-field inflation models induced by the F-term 



First we consider the case inflation is driven by a F-term according to a small-field scenario. We denote the fields 
in the inflaton sector by / collectively although this sector usually contains two or more fields. We assume a minimal 
Kahler potential Ki and a superpotential Wi that we do not specify in details. In the following, we define the F-term 
scalar potential of each sector by 



311^, 



(CI) 



where DjWi = djWi + djKWi. Our fundamental assumption is that the inflaton and modulus sectors are separated. 
Technically, this means K = Kj + Kg and W = Wj + Wg- Then, the total potential reads 



.tot - - Li-i-i ^\DsW\^ 
where the interaction potential is given by 



3\W\'] = el^lV/^)(/) + e^'V^iS) + V^\l, S) , 
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-\S\^\Wi\'^ - ZWjWg - ZWiWl 
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"^here 



Note that the energy scales of Vi and Vg are given by V/(/) ~ 3H]Mp^ and Vg{S) 
Hj is the Hubble parameter during inflation. Let us also remark that the structure of Eq. (|C2|) is in fact typical 
of supergravity [80|: although the two sectors are separated, there is an interaction between them because they 
"communicate" through gravity. In the Minimal Supersymmetric Standard Model (MSSM), this is a possible way to 
transmit the breaking of supersymmetry to the observable sector ^fl| . The same mechanism has also been used to 
argue that the quintessence and inflaton fields must interact (sH or to compute the interaction of dark energy with 
the observable sector, see Ref. [82|. 

Since we are assuming a small-field inflation here one can easily see that jU/W/l > |W/|, and Vi « |-D/W7p « 
3HfMpy In this case the modulus acquires a mass ~ Vj/Mp^ « 3Hf during inflation. After inflation, |/| takes a value 
^ Mpi and oscillates around the potential minimum / = ly^in with Vf(/min) = 0, so that Wi and DjWi are of the 
same order of magnitude on average. If the field oscillation is dominated by a quadratic potential, we find the time 
average of the potential energy Vi{I) = ptot{t)/2 — 3_ff^(i)Mp[/2, so that we can parameterize DjWi = ciHMpi, 
Wi = C2i?Mp[ with ci and C2 parameters of order of unity. If we take Ki = |/ — /minP, the effective potential for S 
averaged over an inflaton oscillation period reads 
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Here /i and /2 are numerical coefficients defined by /i — 2/3(c2 -|- C2)e^^^ and /2 — 2c2e^~^, respectively. Thus in 
the early field oscillation stage when H is larger than the modulus has a mass larger than H and its minimum 

is located at 



f2m3/2MpiH 



(3/2 + 4) - hm3/2H 



(C5) 



As H decreases, it eventually settles down to the absolute minimum S'min — (v^ — ^)Mpi with a mass •\/2m3/2 
determined by Vg{S). Hence this model can be regarded as an example of models where the modulus acquires a mass 
larger than the Hubble parameter during both inflation and the subsequent fleld oscillation regimes with a shift of 
the minimum of order of Mpi typically. 
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2. D-term hybrid inflation 



Next we consider the case in which inflation is induced by the D-term potential. As the simplest realization we 
consider a hybrid inflation model originally proposed in Refs. [H, [sj. Although this model is not observationally 
viable any more due to the largeness of the cosmic string tension produced after inflation, a number of remedies 
have been proposed e.g. in Ref. [s^. In order to avoid inessential complexity, we stick to the original model here 
to adopt the superpotential with three chiral superfields W/ = XQ^i+cj)^ together with Ws — Mg{f3 + S) as in the 
previous subsection. (j)± has a U(l) gauge charge ±1, while Q is neutral. Assuming the minimal Kahler potential. 
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+ 10- P + jS'P, the scalar potential is given by Vtc 



V^"\ with 
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f + \Q4>-? + |0+</'-P + 3|Q0+0_|2) + (|0+|2 + |0_|2 + |^|2 + |Q|2) |AQ0+0_ + Af|(/3 + ^)| 



-3M|/32 + Ml _ 2M|/3(S' + S^) 



(C6) 



where the last term is the D-term contribution with ^ > being the Fayet-IUiopoulos (FI) term. For \Q\> gV^/X = 
Qcri, the potential minimum with respect to (j)± is found at 0± = and inflation is induced by the D-term energy 
density jl. Without Wg the F-term potential vanishes in this regime and the motion of the inflaton is governed 
by a potential generated by quantum corrections. Then terms of the potential relevant to S are given by 



(S, g, 0± = 0) = Mle\^l''+lQ\'' + igp) 1/3 + 5p - 3/3^ - 2/3 (5 + ^t) _ \s\^ + i] 



so that the modulus acquires an extra mass-squared. 



Ml 



J, ,r2 "^3/2 > 



(C7) 
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which is smaller than the original one as long as \Q\ is smaller than Mpi. The shift of the potential minimum of 5* is 
of order of IQp/Mpi. Thus in this model the modulus mass remains much smaller than the Hubble parameter during 
inflation, and the shift of the minimum may also be much smaller than Mpi depending on the value of Q. Note, 
however, that the field configuration of S is determined by long-wave quantum fiuctuations in this case. 

As Q becomes smaller than Qcri the instability occurs with respect to (j>- and infiation is terminated. In this regime 
the potential for S acquires an additional term, SVs{S) = e''^' X'^\Q(j)-\'^ . This term induces a correction to the 
modulus mass term 



SeS — 



Ml, 



9 



(C9) 



where we have inserted the minimum of for Q < Qai- The last expression takes its maximum at \Q\'^ = |QcriP/2 
and is given by Sm'i.^ff — g^£^'^ /A. This means that after infiation the effective mass of the Polonyi field rises significantly 
to the level comparable to the Hubble parameter at that time. When this induced mass term is operative, the potential 

minimum for \S\ takes a value O (^m^.^/ H^^ Mpi which departs from the eventual minimum by Mpj. 



3. Hybrid inflation with both F- and D-term contributions 



Next we consider a variant of D-term hybrid inflation model with a non- vanishing F-term potential as an example 
of models where the modulus mass remains much smaller than the Hubble parameter during infiation but it acquires 
a lar ge c orrection of order H just after infiation. Specifically we consider sneutrino hybrid inflation model proposed in 
Ref. I86| with the minimal Kahler potential. The effects of non-minimal Kahler potential are discussed, for example, 
in Ref. |87| for the sneutrino inflation and in Ref. (ssj for D-term infiation. 

The infiaton sector of the model consists of three species of chiral superfields, containing (s)neutrinos (i = 1,2,3), 
and (f)± with U(l) gauge charge ±1 as in the above model. The relevant part of the superpotential reads 

Wi = ^iVfiVf + ^M.N^N^ , (CIO) 
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besides the interaction between iVf and lepton and Higgs fields. One can always choose a basis for so that their 
mass matrix is diagonal, and we take to be Majorana mass eigenstate fields with real mass M^. We assume, 
without loss of generality, the inflaton sneutrino is the lightest heavy sneutrino Mi <C M2 , M3 and we are interested 
in the lower range of the preferred values of heavy neutrino masses Adi — ~ lO^^GeV. D-term contribution to 
the scalar potential with non- vanishing FI term ^ > is given by Vd 0_) in Eq. (|C6p as before. 

The large amplitude of inflaton sneutrino gives a large effective mass to the slepton and the Higgs field Hu and 
they stay at the origin and do not affect the inflation dynamics. 04- has a positive mass during and after inflation 
and it stays at the origin all the time. For the discussion of inflationary dynamics, therefore, we discuss the evolution 
of iVi and . The tree-level potential for the inflaton sector is then given by 



+ = e^^ 



(Cll) 



We expect A 3> M/Mpi, and therefore find 



for 



> g^^Mpj/A^ = A'^j.j during inflation. (/>- destabilizes 



for iVi < A^cri to reach its minimum at 0_ — after inflation. We find A'cri ^ 2Mp\ for cosmologically relevant 
parameter values [86| . 

The correction to the Polonyi mass during inflation can be found from the F-term potential with (f)± = 0. The 
interaction term reads 
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For iVi ^ Mp\ we find the last term has the largest contribution to the effective modulus mass, 0{M). This correction, 
however, is still much smaller than the Hubble parameter during inflation thanks to the assumption that inflation is 
driven by the D-term potential energy which ensures that 



(C13) 



As Ni gets smaller than Ncn, a phase transition occurs to make nonvanishing and inflation is terminated. Its 
minimum is located at 



(C14) 



Wc therefore find 



where the equality holds at 
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/2 in the last inequality. Thus one can see that in the early field oscillation 



regime after inflation, the modulus acquires an effective mass of order of the Hubble parameter during inflation. 
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